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ABSTRACT 

We present a theoretical interpretation of the broad silicon and iron ultraviolet ab- 
sorption features observed with the Hubble Space Telescope in the spectrum of the 
Schweizer-Middleditch star behind the remnant of Supernova 1006. These features are 
caused by supernova ejecta in SN1006. We propose that the redshifted Si ll 1260 A 
feature consists of both unshocked and shocked Si ll. The sharp red edge of the line 
at 7070 km s~^ indicates the position of the reverse shock, while its Gaussian blue edge 
reveals shocked Si with a mean velocity of 5050 km s~^ and a dispersion of 1240 km s~^, 
implying a reverse shock velocity of 2860 km s~^. The measured velocities satisfy the 
energy jump condition for a strong shock, provided that all the shock energy goes into 
ions, with little or no collisionless heating of electrons. The line profiles of the Si ill 
and Si IV absorption features indicate that they arise mostly from shocked Si. The total 

^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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mass of shocked and unshocked Si inferred from the Si ll, Si ill and Si IV profiles is 
Msi = 0.25 ± 0.01 Mq on the assumption of spherical symmetry. Unshocked Si extends 
upwards from 5600 km s~^. Although there appears to be some Fe mixed with the Si 
at lower velocities ^ 7070 km s~^, the absence of Fe ll absorption with the same profile 
as the shocked Si ll suggests little Fe mixed with Si at higher (before being shocked) 
velocities. The column density of shocked Si ll is close to that expected for Si ll un- 
dergoing steady state coUisional ionization behind the reverse shock, provided that the 
electron to Si II ratio is low, from which we infer that most of the shocked Si is likely 
to be of a fairly high degree of purity, unmixed with other elements. We propose that 
the ambient interstellar density on the far side of SN1006 is anomalously low compared 
to the density around the rest of the remnant. This would simultaneously explain the 
high velocity of the redshifted Si absorption, the absence of blueshifted Si absorption, 
and the low density of the absorbing Si compared to the high Si density required to 
produce the observed Si x-ray line emission. We have reanalyzed the Fe ll absorption 
lines, concluding that the earlier evidence for high velocity blueshifted Fe II extending 
to ~ —8000 km s~^ is not compelling. We interpret the blue edge on the Fe ll profiles at 
—4200 km s~^ as the position of the reverse shock on the near side of SN1006. The mass 
of Fe 11 inferred from the red edge of the Fe ll profile is Mpeii = 0.029 ± 0.004 Mq up to 
7070 km s~^, if spherical symmetry is assumed. The low ionization state of unshocked 
Si inferred from our analysis of the silicon features, Si Il/Si = 0.92 ± 0.07, suggests a 
correspondingly low ionization state of unshocked iron, with Fe il/Fe = 0.66]!Iq;22- If this 
is correct, then the total mass of Fe up to 7070 km s~^ is Mpe = 0.044lo:oi3 Mq with 
a 3(7 upper limit of Mpe < 0.16 Mq. Such a low ionization state and mass of iron is 
consistent with the recent observation of Fe III 1123 A with HUT, indicating Fe IIl/Fe II 
= 1.1 ± 0.9, but conflicts with the expected presence of several tenths of a solar mass of 
iron in this suspected Type la remnant. However, the inference from the present HST 
data is too indirect, and the HUT data are too noisy, to rule out a large mass of iron. 
Re-observation of the Fe III 1123 A line at higher signal to noise ratio with FUSE will 
be important in determining the degree of ionization and hence mass of iron in SN1006. 

Subject headings: ISM: individual (SN1006) — ISM: supernova remnants — ultraviolet: 
ISM 
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1. Introduction 

In two previous papers, we (Wu et al. 1993, 1996, 
hereafter WCFHS93, WCHFLS96) described HST 
FOS observations of the UV spectrum of the SM 
star (Schweizer & Middleditch 1980) which lies be- 
hind and close to the projected center of the remnant 
of SN1006. In the present paper we offer a theoret- 
ical interpretation of the broad sihcon and iron UV 
absorption features observed with HST. These fea- 
tures are almost certainly caused by supernova ejecta 
in SN1006, as originally proposed by Wu et al. (1983), 
who first observed the features with lUE. 

Detailed theoretical analysis of the Fe II features 
observed with lUE has been presented previously by 
Hamilton & Fesen (1988, hereafter HF88). The main 
purpose of that paper was to try to explain the ap- 
parent conflict between the low « 0.015 Mq mass of 
Fe II inferred from the lUE observations of SN1006 
(Fesen et al. 1988) with the expected presence of sev- 
eral tenths of a solar mass of iron (Hoflich & Khokhlov 
1996) in this suspected Type la remnant (Minkowski 
1966; Schaefer 1996). HF88 demonstrated that ambi- 
ent UV starlight and UV and x-ray emission from 
reverse-shocked ejecta could photoionize unshocked 
iron mainly to Fe III, Fe IV, and Fe V, resolving the 
conflict. 

Recently Blair, Long & Raymond (1996) used 
the Hopkins Ultraviolet Telescope (HUT) to measure 
Fe III 1123 A absorption in the spectrum of the SM 
star. They found Fe III/Fe II = 1.1 ± 0.9, which 
neither confirms, nor excludes, the ratio Fe III/Fe II 
= 2.6 predicted by HF88. 

The HST spectra, particularly the silicon features, 
prove to be a rich source of information beyond the 
reach of lUE's capabilities. In the first half of this 
paper, Section ||, we a nalyze the Si absorption fea- 
tures. We find (j j2.l| ) that the profile of the red- 
shifted Si II 1260 A feature, with its sharp red edge 
and Gaussian blue edge, appears to be attributable to 
the presence of both unshocked and shocked silicon. 
We then develop a chain of inferences, first about the 
reverse shock (§^|^) and collisionless heating (^^T^) 



then about the column density and ma ss (§§2.4, 2.7), 
purity (§ 2.5), an d ionization state (j |2.6| ) of the silicon. 
We argue (§2^) that the ambient interstellar density 
on the far of SN1006 is anomalously low compared 
to density around the rest of the remnant, which ex- 
plains the high velocity of the redshifted Si, the ab- 
sence of corresponding blueshifted Si (§ 2.9 ) , and some 



other observational puzzles. 

In the second half of the paper. Section ||, we dis- 
cuss the broad Fe II absorption features. WCFHS93 
reported blueshifted Fe II absorption up to —8000 kms" 
Finding the presence of such high velocity blueshifted 
absorption difficult to understand in the light of other 
observational evidence, we detail a reanalysis of the 
Fe II features in 
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We conclude (§3.2) that the ev- 
idence for high velocity blueshifted absorption is not 
compelling, and we propose ( |3.3| ) that the sharp blue 
edge on the Fe II profiles at — 4200kms^^ represents 
the position of the reverse shock on the near side. In 
the remainder of the Section, §§3.4-3.6, we address 



the issue of the ionization state and mass of the iron. 

We attempt in this paper to construct a consistent 
theoretical picture, but there remain some discrepan- 
cies, and we highlight these in Section ^. Section || 
summarizes the conclusions. 

2. Silicon 

There are three possibilities for the origin of the 
broad Si absorption features in the spectrum of the 
SM star, if it is accepted that these features arise from 
the remnant of SN1006. 

The first possibility is that the absorption arises 
from cool, dense, fast moving knots of ejecta, as sug- 
gested by Fesen et al. (1988) and further discussed by 
Fesen & Hamilton (1988). However, there has been no 
significant change in the features over 12 years, from 
the original detection of the features with lUE in 1982 
(Wu et al. 1983), through their re-observation with 
lUE in 1986 and 1988 (Fesen & Hamilton 1988), up to 
the FOS observation with HST in 1994 (WCHFLS96). 
Furthermore, the relative strengths of the redshifted 
Si II 1260 A, 1304 A, and 1527 A features are approxi- 
mately proportional to their oscillator strengths times 
wavelengths, indicating that the lines are not satu- 
rated. This constancy in time and lack of saturation 
argues against the absorption features being caused 
by small, dense knots. We do not consider this hy- 
pothesis further in this paper. 

A second possibility is that the Si absorption is 
from shocked ejecta in which the coUisional ioniza- 
tion timescale is so long that the observed low ion- 
ization species Si II, Si III, and Si IV can survive. 
At first sight, the high ^ 5000 kms~^ velocity of the 
observed absorption makes this possibility seem un- 
likely, because shocked ejecta should be moving no 
faster than the velocity of gas behind the interstellar 
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shock, which in the NW sector of the remnant can be 
inferred from the 2310 kms^^ FWHM of the Balmer 
broad hne emission to be 1800-2400 kms~^ (this is 
3/4 of the shock velocity), depending on the extent 
to which electron-ion equilibration takes place (Kir- 
shner, Winkler & Chevalier 1987; Long, Blair & van 
den Bergh 1988; Smith et al. 1991; Raymond, Blair & 
Long 1995). However, below we wiU conclude that it 
is likely that much, in fact most, of the Si absorption 
is from shocked ejecta, and that the ISM surrounding 
SN1006 may be quite inhomogeneous. 

A third possibility is that the Si absorption arises 
from unshocked supernova ejecta freely expanding in 
SN1006, which is consistent with the high velocity of 
the absorption. The low ionization state of the Si, 
predominantly Si II, with some Si III and Si IV, is at 
least qualitatively consistent with the expectations of 
models in which the unshocked ejecta are photoion- 
ized by ambient UV starlight and by UV and x-ray 
emission from shocked ejecta (HF88). Neutral Si is 
neither observed, e.g. at Si I 1845 A, nor expected, 
since it should be quickly 20 yr) photoionized by 
ambient UV starlight. Recombination is negligible at 
the low densities here. At the outset therefore, this 
possibility seems most likely, and we pursue the idea 
further in the next subsection. 

Fesen et al. (1988) pointed out that the redshifted 
Si II 1260 A feature (at - 1280 A) in the lUE data 
appeared to be somewhat too strong compared to the 
weaker redshifted Si II 1527 A feature. The discrep- 
ancy appears to be confirmed by the HST observa- 
tions (WCHFLS96). Fesen et al. proposed that some 
of the Si II 1260 A feature may come from S II 1260, 
1254, 1251 A redshifted by « 800 kms^^ relative to 
the Si II, a possibility also addressed by WCHFLS96. 
In the present paper we regard the possibility of any 
significant contribution from S II as unlikely, notwith- 
standing the discrepancy between the Si II profiles. 
The oscillator strengths of the S II 1260, 1254, 1251 A 
lines are / = 0.01624, 0.01088, 0.005453, whose com- 
bined strength is only 1/30 of the oscillator strength 
/ = 1.007 of the Si II 1260 A line (Morton 1991). 
In models of Type la supernovae, such as SN1006 is 
believed to be, silicon and sulfur occur typically in 
the same region of space, with a relative abundance 
of Si : S w 2 : 1 (e.g. Nomoto, Thielemann & Yokoi 
1984). Thus S II might be expected to contribute only 
~ 1/60 of the optical depth of Si II in the 1260 A fea- 
ture, assuming a similar ionization state of Si and S. 
In the remainder of this paper we ignore any possible 



contribution of S II to the Si II 1260 A feature. In this 
we follow Wu et al.'s (1983) original identification of 
the 1280 A absorption feature as redshifted Si. 

2.1. The redshifted Si II 1260 A feature 

Given that SN1006 shows a well-developed inter- 
stellar blast wave in both radio (Reynolds & Gilmore 
1986, 1993) and x-rays (Koyama et al. 1995; Will- 
ingale et al. 1996), it is inevitable that a reverse shock 
must be propagating into any unshocked ejecta. 

If the observed UV Si absorption is from unshocked 
ejecta, then there should be a sharp cutoff in the 
line profile at the expansion velocity of the reverse 
shock, because the shock should 'instantaneously' de- 
celerate the ejecta to lower bulk velocity. In fact the 
Si II 1260 A feature does show a steep red edge at 
7070 ± 50kms^^, albeit with a possible tail to higher 
velocities. Tentatively, we take the presence of the 
steep red edge as evidence that at least some of the 
Si II is unshocked. 

Once shocked, how long will Si II last before being 
collisionally ionized to higher levels? The ionization 
timescale of Si II entering the reverse shock can be 
inferred from the optical depth of the Si II 1260 A ab- 
sorption just inside its steep red edge at 7070 kms~^, 
which putatively represents the position of the reverse 
shock. In freely expanding ejecta where radius equals 
velocity times age, r = vt, the column density per 
unit velocity dN/dv of any species is equal to the 
density n — dN/dr times age t. For freely expanding 
ejecta, the optical depth r in a line of wavelength A 
and oscillator strength / is then proportional to nt: 

ne^ dN ne^ 
T = A-l- = f>'n-t . (1) 

TOgC dv TUeC 

The optical depth of the Si II 1260.4221 A line (/ = 
1.007, Morton 1991) just inside its steep red edge is 
T w 1. This implies, from equation (|l|), a preshock 
Si II density times age of ?^g[n '^^ — 3.0 x lO^cm^'^s. 
The postshock Si II density times age would then 
be 4 times higher for a strong shock, nsiii^ = 1.2 x 
10^cm~'^s. At a coUisional ionization rate for Si II 
of (crw)siii = 6 X lO'^cm^s^i (Lennon et al. 1988; 
see subsection below for further discussion of 
this rate), the ratio of the Si II ionization timescale 
igiii = {ne{<yv)siu)^^ to the age of the remnant, t, is 
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Table 1 

Parameters measured from redshifted Si II 1260 A feature 



Parameter 


Value 




Expansion velocity into reverse shock 


7070 ±50 km s" 


-1 


Mean velocity of shocked Si II 


5050 ±60 km s- 


-1 


Dispersion of shocked Si II 


1240 ±40 km s" 


-1 


Reverse shock velocity 


2860± 100kms- 


-1 


Lower edge of unshocked Si II 


5600± lOOkms" 


-1 


Preshock density of Si II 


5.4 ±0.7 X 10^5 cm- 


-3 


Column density of shocked Si II 


9.0 ±0.3 X 10^^ cm' 


-2 


Column density of unshocked Si II 


1.5 ±0.2 X IQi-^cm- 


-2 


Column density of all Si II 


10.5 ±0.1 X IQi'^cm- 


-2 


Mass of shocked Si II 


0.127 ± 0.006 M 





Mass of unshocked Si II 


0.017 ± 0.002 M 






Note. — Masses assume spherical symmetry. For simplicity, 
velocities and masses have not been adjusted for the small offset 
of the SM star from the projected center of the remnant. 



(the following estimate is revised below, equation g]) 

tsai _ 1 _ '^Siii 1 _ ^ ^ "Siii 

t net{av)siu Tie nsiiii((Tw)siii ' 

(2) 

Since the ratio ne/n^iii of electron to Si II density in 
the postshock gas should be greater than but of order 
unity, this estimate (||) indicates that the coUisional 
timescale of Si II is of the order of the age of the 
remnant. It follows that shocked Si II is likely also to 
contribute to the observed Si II absorption. 

Shocked Si II will be decelerated by the reverse 
shock to lower velocities than the freely expanding 
unshocked ejecta. Shocked Si II should have a broad 
thermal profile, unlike the unshocked Si II. Examin- 
ing the redshifted Si II 1260 A profile, we see that 
the blue edge extends down to about ±2500 kms~^, 
with a shape which looks Gaussian. Fitting the blue 
edge to a Gaussian, we find a best fit to a Gaussian 
centered at 5050 kms~^, with a dispersion (standard 
deviation) of cr = 1240 kms^^. This fit is shown in 
Figure |l|. Having started from the point of view that 
the Si II was likely to be unshocked, we were sur- 
prised to see that, according to the fit, it is shocked 
Si II which causes most of the absorption, although an 



appreciable quantity of unshocked Si is also present, 
at velocities extending upwards from 5600 km s~^. 
The slight tail of Si II absorption to high velocities 
> 7070 km s~^ is naturally produced by the tail of 
the Gaussian profile of the shocked Si II. 

The estimate (||) of the coUisional ionization timescale 
of Si II presumed that all the Si II 1260 A absorption 
was from unshocked Si II, whereas the picture now is 
that only some of the absorption is from unshocked 
Si II. According to the fit in Figure |^, the optical 
depth of unshocked Si II at the reverse shock front 
is T = 0.56 ± 0.07, a little over half that adopted in 
estimate (^, so a revised estimate of the ionization 
timescale of Si II is not quite double that of the orig- 
inal estimate (^): 

t He 

Evidently the conclusion remains that the ionization 
timescale of Si II is comparable to the age of the rem- 
nant. 

2.2. Shock jump conditions 

The fitted profile of the Si II 1260 A feature in 
Figure H includes both unshocked and shocked com- 
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ponents. The consistency of the fitted parameters 
can be checked against the jump conditions for a 
strong shock. The shock jump conditions predict that 
the three-dimensional velocity dispersion of the 

ions should be related to the deceleration Av of the 
shocked gas by energy conservation 



= Av 



(4) 



provided that all the shock energy goes into ions. The 
observed dispersion is 

= X (1240±40kms"^) = 2140±70kms-^ 

(5) 

while the observed deceleration is 

Av = (7070±50kms~^) - (5050±60kms"^) 
= 2020±80kms"i . (6) 

These agree remarkably well, encouraging us to be- 
lieve that this interpretation is along the right lines. 
The reverse shock velocity, v^, corresponding to the 
observed dispersion is 



= (16/3)^/^ = 2860 ± 100 kms" 



(7) 



We prefer to infer the shock velocity from the ob- 
served dispersion rather than from the observed de- 
celeration Av, since the latter may underestimate the 
true deceleration, if, as is likely, the shocked Si is 
moving on average slightly faster than the immedi- 
ate postshock gas (see below). 

The predicted equality between the decelera- 
tion and ion dispersion holds provided that all the 
shock energy goes into ions, and that the bulk veloc- 
ity and dispersion of the ions in the shocked gas are 
equal to their postshock values. Since the shocked 
Si II can last for a time comparable to the age of the 
remnant (equation |^]), it cannot be assumed auto- 
matically that the bulk velocity and dispersion of the 
observed Si II ions are necessarily equal to those im- 
mediately behind the reverse shock front. We discuss 
first the issue of the bulk velocity, then the disper- 
sion, and fin ally the question of coUisionless heating 
in subsection 2.3. 



Consider first the bulk velocity of the shocked 
ions. In realistic hydrodynamic models, the veloc- 
ity of shocked gas increases outward from the re- 
verse shock. Indeed, the fact that the observed dis- 
persion 3^/^cr is larger than the observed decelera- 
tion Av by 120 ± 130 kms^^ (the uncertainty here 



takes into account the correlation between the un- 
certainties in a and Av) is consistent with the no- 
tion that the shocked Si II is moving on average 
120 ± 130kms~^ faster than the immediate post- 
shock gas. This modest velocity is consistent with ex- 
pectations from one-dimensional hydrodynamic sim- 
ulations appropriate to SN1006 (see HF88, Fig. 2), 
according to which shocked ejecta remain relatively 
close to the reverse shock front. However, the de- 
celeration of ejecta is generally Rayleigh- Taylor un- 
stable (e.g. Chevalier, Blondin & Emmering 1992), 
which instabilities could have caused shocked Si II 
to appear at velocities many hundred kms~^ faster 
than the immediate postshock gas. Since the obser- 
vations do not show this, it suggests, though by no 
means proves, either that Rayleigh- Taylor instabili- 
ties are not very important, or perhaps that the line 
of sight through SN1006 to the SM star happens to 
lie between Rayleigh- Taylor plumes. 

What about the ion dispersion? If there were a 
range of ion dispersions in the shocked gas, then the 
line profile would tend to be more peaked and have 
broader wings than a simple Gaussian. The observed 
profile of the blue edge of the redshifted Si II 1260 A 
feature is consistent with a Gaussian, which suggests, 
again weakly, that the ion dispersion does not vary 
by a large factor over the bulk of the shocked Si II. 

While the observations agree well with the simplest 
possible interpretation, it is certainly possible to ar- 
range situations in which a combination of Rayleigh- 
Taylor instabilities, spatially varying ion dispersion, 
and coUisionless heating conspire to produce fortu- 
itous agreement of the observations with the jump 
condition (||). 

2.3. CoUisionless heating 

The shock jump condition (Q) is valid provided that 
all the shock energy goes into the ions, rather than 
into electrons, magnetic fields, or relativistic particles. 
The timescale for equilibration by Coulomb collisions 
between electrons and ions or between ions and ions 
is much longer than the age of SN1006. However, 
coUisionless processes in the shock may also transfer 
energy, and the extent to which these may accom- 
plish equilibration remains uncertain (see Laming et 
al. 1996 for a recent review). The prevailing weight 
of evidence favors little if any coUisionless heating of 
electrons in the fast shocks in SN1006. Raymond et 
al. (1995) find similar velocity widths in emission lines 
of H I (Ly /3), He II, C IV, N V, and O VI observed by 
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HUT from the interstellar shock along the NW sector 
of SN1006. They conclude that there has been lit- 
tle equilibration between ions, though this does not 
rule out substantial electron heating. From the same 
data, Laming et al. (1996) argue that the ratio of C IV 
(which is excited mainly by protons) to He II (which is 
excited mainly by electrons) is again consistent with 
little or no electron-ion equilibration. Koyama et al. 
(1995) present spectral and imaging evidence from 
ASCA that the high energy component of the x-ray 
spectrum of SN1006 is nonthermal synchrotron ra- 
diation, obviating the need for coUisionlcss electron 
heating (see also Reynolds 1996). 

The results reported here, equations (|^) and (|^), 
tend to support the conclusion that virtually all the 
shock energy is deposited into the ions. If some of the 
shock energy were chaneled into electrons, magnetic 
fields, or relativistic particles, then the ion dispersion 
would be lower than predicted by the observed decel- 
eration, whereas the opposite is observed — the ion 
dispersion is slightly higher. 

2.4. Column density and mass of Si II 

The numbers given here are summarized in Ta- 
ble 1^. Quoted uncertainties in column densities and 
masses here and throughout this paper ignore uncer- 
tainties in oscillator strengths. The uncertainties in 
the oscillator strengths of Si II 1260 A is 17%, and of 
Si III 1206 A and Si IV 1394, 1403 A are 10% (Morton 
1991). 

The column density of shocked and unshocked Si II 
follows from integrating over the line profiles shown 
in Figure |^. The column density Ngf^^^ of shocked Si II 
is 

AT^h ^ ^!^(2^)i/2^^p ^ 9.0±0.3xlOi4cm-2 (8) 

where tq = 0.98 ±0.02 is the optical depth at line cen- 
ter of the fitted Gaussian profile of the shocked Si II, 
and the factor (27r)^''^cr comes from integrating over 
the Gaussian profile. The profile of unshocked Si II is 
the residual after the shocked Si II is subtracted from 
the total, and we measure its column density N^^f^ 
by integrating the unshocked profile over the velocity 
range 5600-7200 kms^i; 

N^Rf = / T""^'^dt; = 1.5±0.2xl0i*cm-2 (9) 

TTe^fXJ 

where the uncertainty is largely from uncertainty in 
the fit to the shocked Si II. Integrating the full Si II 



profile over 1000-8200 kms-i yields the total column 
density of Si II 

iVsiii = 10.5 ± 0.1 X 10"cm-2 (10) 

where the uncertainty is from photon counts, and is 
smaller than the uncertainties in either of the shocked 
or unshocked Si II column densities individually (be- 
cause there is some uncertainty in allocating the total 
column density between the shocked and unshocked 
components). 

The corresponding masses of unshocked and shocked 
Si II can also be inferred, if it is assumed that the 
Si was origin ally ejected spherically symmetrically. 
In subsection ^.8| we will argue that the absence of 
blueshifted Si II absorption can be explained if the 
reverse shock has passed entirely through the Si on 
the near side of SN1006, and Si has been collision- 
ally ionized to high ion stages. Thus the absence of 
blueshifted Si II need not conflict with spherical sym- 
metry of Si ejected in the supernova explosion. If 
the shocked ejecta are taken to lie in a thin shell at 
a free expansion radius of v — 7200 kms^^ (slightly 
outside the position of reverse shock — cf. the argu- 
ment in paragraph 3 of subsection |]^), then the mass 
of shocked Si II is 

M^lii = 47rmsiii(wt)^iVsni 0.127 ± 0.006 Mq (11) 

where the uncertainty includes only the uncertainty 
in the column density of shocked Si II. This value 
should perhaps be modified to M^j^jj = 0.13±0.01 Mq 
to allow for uncertainty in the radial position of the 
shocked Si II. The mass of unshocked Si II is an inte- 
gral over the unshocked line profile 

M^Rf = ATTmsiut^J v^dN^Rf = 0.017 ± 0.002 Mq 

(12) 

where the uncertainty is largely from uncertainty in 
the fit to the shocked Si II. 

2.5. Purity of shocked Si 

We show below that the observed column density 
of Si II is close to the column density which would be 
predicted under the simple assumption of steady state 
collisional ionization downstream of the shock (note 
that recombination is negligible). Now the assump- 
tion of steady state ionization is surely false, since the 
timescale to ionize Si II is comparable to the age of the 
remnant, equation (0). Below we will argue that the 
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effect of non-steady state is generally such as to reduce 
the column density below the steady state value. The 
constraint that the observed column density should 
be less than or comparable to the steady state value 
then leads to an upper limit on the electron to ion 
ratio, equation (14). The fact that this upper hmit 



is not much greater than unity leads to the interest- 
ing conclusion that the bulk of the observed shocked 
Si II must be fairly pure, since admixtures of other el- 
ements would increase the electron to ion ratio above 
the limit. 

If the postshock number density of Si II ions is 
nsiii (which is 4 times the preshock number density), 
then at shock velocity Vg the number of Si II ions 
entering the reverse shock per unit area and time is 
nsiuVs/4:. The Si II ions are collisionally ionized by 
electrons in the shocked gas at a rate ?T.e(o'w)siii ion- 
izations per unit time. For Vg — 2860 kms~^ and 



{crv)siu = 6.1 X 10" 



^cm'^s 



it follows that the col- 



umn density A^g^J^j'''*^ of shocked Si II in steady state 
should be 

^steady ^ nsiUVg ^ ^ ^ x lOiScm^^ . (13) 

If, as will be argued below, the actual (non-steady 
state) column density (||) of shocked Si II is less than 
or comparable to the steady state value (p^, then the 
mean ratio of electron to Si II density in the shocked 
gas satisfies 

^ 1.3 . (14) 

nsiu 

The ratio must of course also satisiy ne/nsm > 1, 
since each Si II ion itself donates one electron. Such a 
low value ( p^ of the electron to Si II ratio in the 
shocked Si II would indicate that the bulk of the 
shocked Si II must be of a rather high degree of purity, 
since the presence of significant quantities of other el- 
ements or higher ionization states of Si would increase 
the number of electrons per Si II ion above the limit 
(14). Indeed, even if the Si II entering the shock were 
initially pure, ionization to Si III and higher states 
would release additional electrons. In steady state, 
the mean electron to Si II ratio experienced by Si II 
during its ionization from an initially pure Si II state 
is ne/nsiii = 1.5, already slightly larger than the limit 



The limit ( [T^ ) on the electron to Si II ratio is so low 
as to make it difficult to include even modest quanti- 
ties of other elements with the silicon. This may be 
a problem. While Si is generally the most abundant 



element in Si-rich material produced by explosive nu- 
cleosynthesis, other elements, notably sulfur, usually 
accompany the silicon. For example, the deflagrated 
white dwarf model W7 of Nomoto et al. (1984) con- 
tains 0.16 Mq of Si mostly in a layer which is about 
60% Si, 30% S by mass. At this elemental abun- 
dance, and assuming similar ionization states for all 
elements, the expected mean electron to Si II ratio 
in steady state would be rie/nsiii = 1.5/0.6 = 2.5, 
almost twice the limit given by equation (p^). Be- 
cause of this potential difficulty, we discuss carefully 
below how robust is the constraint (|lj) . First we ad- 
dress the accuracy of the predicted value (13) of the 
steady state column density, and then we discuss the 
non-steady state case. 

The electron to Si II ratio could be increased if 
the predicted steady state column density ( [l3| ) were 
increased, either by increasing the shock velocity Vs, 
or by reducing the coUisional ionization rate (cru)siii. 
We consider the former first, then the latter. The 
shock velocity Vg = 2860 km s~^ is inferred from the 
observed ion dispersion in the Si II 1260 A line, equa- 
tion (0). This should give a fair estimate of the mean 
shock velocity of the observed shocked Si II, except 
for a small correction from the fact that the ion dis- 
persion (temperature) in the past would have been 
higher because of adiabatic expansion of the rem- 
nant. Moffet, Goss & Reynolds (1993) find that the 
global radius R of the radio remnant is currently in- 
creasing with time according to i? oc iO-48±o.i3 j£ 
the ambient ISM is assumed to be uniform, this in- 
dicates that the pressure in the remnant is varying 
with time as P cx {R/t)"^ oc t-i o4±o.26^ hence that the 
temperature is varying in Lagrangian gas elements as 
T oc P^/^ oc i-0-42±o.io^ hence that the ion dispersion 
is varying as cr oc T^/^ oc i-0-2i±o.05^ rather weak 
function of time. If one supposes that the observed 
Si II was shocked on average when SN1006 was say 
half its current age, then the dispersion, hence the 
reverse shock velocity, could have been ~ 20% higher 
than at present. The steady state column density 
( |l3| ) would then be 20% higher, and the constraint 
on the electron to Si II ratio (^^ would be relaxed 
slightly to rie/nsiii <; 1.5. 

The collisional ionization rate (crw)siii — 6.1 x 
10~^cm'^s~^ used in equation ([l3| ) comes from inte- 
grating the cross sections of Lennon et al. (1988) over 
a Maxwellian distribution of electrons at a tempera- 
ture of 83 eV. The quoted error in the cross-sections is 
60%, the large uncertainty arising from the fact that 
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the cross-sections for Si II are derived from isoelec- 
tronic scaling rather than from real data. Reducing 
the ionization rate by 0.2 dex would relax the con- 
straint (|l^) to ne/nsiii < 2.0. 

The temperature of 83 eV used in the coUisional 
ionization rate above is the temperature reached by 
electrons as a result of Coulomb collisions with Si II 
ions over the coUisional ionization timescale of Si II. 
The assumption here that there is no coUisionless 
heating of electrons in the shock is in accordance 
with the arguments given in subsection 2.3. Actu- 



ally, the ionization rate of Si II by electron impact, 
as derived from the cross-sections of Lennon et al., 
has a broad maximum at an electron temperature of 
^ 200 eV, and varies over only 5-7 x lO^^cm'^s^^ for 
electron temperatures 40-1000 eV. Thus uncertainty 
in the electron temperature does not lead to much 
uncertainty in the coUisional ionization rate, unless 
there is substantial coUisionless heating of electrons 
to temperatures much higher than 1 keV. 

We have argued against significant coUisionless 
electron heating in subsection 2.3. However, if coUi- 
sionless electron heating to temperatures greater than 
1 keV did occur, it would imply both a higher shock 
velocity, since the observed ion dispersion would un- 
derestimate the shock energy, and a lower coUisional 
ionization rate, both of which act to increase the 
steady state column density (|T^). Thus coUisionless 
electron heating, if it occurs, would allow a larger 
electron to Si II ratio than given by equation . 

We now turn to the argument that in a non-steady 
state situation, as here, the column density of shocked 
Si II is likely to be less than the steady state column 
density (p!^), which leads to the constraint (Q) on 
the mean electron to Si II ratio in the shocked Si II. 
In the first place, simply truncating the Si at some 
point downstream of the shock will give a lower col- 
umn density than the steady state value. Secondly, 
geometric effects tend to reduce the column density 
below the steady state value. That is, the column 
density of Si II is diluted by the squared ratio (r^/r)^ 
of the original radius rg of the gas at the time it was 
shocked to the present radius r (> r^) of this shocked 
gas. Thirdly, if the density profile of shocked Si at the 
present time increases outwards, then the faster ion- 
ization of the denser, earlier shocked, gas reduces its 
column density per interval of ionization time, and 
the net column density is again lower than steady 
state (a more rigorous demonstration of this is given 
in the Appendix) . Conversely, if the density profile of 



shocked Si decreases outwards, then the net column 
density can be higher than steady state, but only if 
the flow is allowed to continue for sufficiently longer 
than a coUisional ionization time. However, according 
to equation (|^) the ionization timescale at the present 
Si II density is comparable to the age of the remnant, 
and the ionization timescale would be longer at lower 
density, so there is not much room for increasing the 
column density this way either. 

So is there any way that the actual column density 
of Si II could be higher than the steady state value? 
Clearly yes, given sufficient freedom with the density 
profile of shocked Si. For example, one possibility is 
that there is a 'hump' in the shocked Si II density 
profile, such that the density increases outward of the 
present position of the reverse shock front, but then 
declines at larger radii. Some tuning is required to 
ensure that the density on both near and far sides 
of the hump is high enough to produce significant 
column density, but not so high as to ionize Si above 
Si II. The higher the column density, the more fine- 
tuning is required. 

We thus conclude that while some violation of the 
limit dl^ ) on the mean electron to Si II ratio in the 
shocked Si II is possible, greater violations, exceeding 
say a factor of two, are less likely. It then follows that 
the bulk of the shocked Si II is likely to be of a fairly 
high degree of purity. In particular, there is unlikely 
to be much iron mixed in with the shocked silicon, 
a conclusion which is consistent with the absence of 
Fe II absorption with the same profile as the shocked 
Si II, as discussed in subsection 3.4, To avoid misun- 



derstanding, this statement refers only to the shocked 
Si II: iron could be mixed with the unshocked Si, and 
indeed the absorption profile of Fe II, Figure || be- 
low, does suggest that there is some Fe mixed with 
unshocked Si. 

2.6. Si III and Si IV line profiles 

Given that shocked Si II apparently persists for a 
time comparable to its ionization time, it is difficult to 
avoid producing an appreciable quantity of Si III and 
Si IV as the result of coUisional ionization of Si II in 
the shocked ejecta. We thus conclude that it is likely 
that most of the observed Si III and Si IV absorption 
arises from shocked ejecta. This is consistent with the 
observed line profiles, as will now be discussed. 

Figures || and || show fits to the redshifted Si III 
1206 A and Si IV 1394, 1403 A features using as tem- 
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plates the shocked and unshocked profiles of Si II 
1260 A shown in Figure 0. Table | gives the fitted 
column densities of shocked and unshocked Si III and 
Si IV, expressed relative to the best fit column density 
of shocked and unshocked Si II given in Table |l|. 

The Si III 1206 A profile appears to be mostly 
shocked. There is some indication of unshocked Si III 
over 5500-7000 kms-i at the 2a level, Table |, as 
suggested by the residual profile after subtraction of 
shocked Si III plotted in Figure ^. The dispersion of 
the fitted Gaussian profile of the Si III, if allowed to 
be a free parameter, is 1290±60kms~^ if unshocked 
Si III is excluded, or 1210 ± 60kms^^ if unshocked 
Si III is admitted, which are in good agreement with 
the 1240 ± 40kms"i dispersion of the Si II. 

The profile of the Si IV 1394, 1403 A feature. Fig- 
ure ^, is consistent with containing no unshocked 
Si IV. If the line center and width of the Gaussian pair 
fitted to the Si IV 1394, 1403 A doublet are allowed to 
be free, then the center remains close to 5050 kms"-'^, 
the same as for Si II 1260 A, but the best fit dispersion 
of the Si IV is 1700± lOOkms^i, which is 4.5ct higher 
than the 1240 kms"^ dispersion of Si II 1260 A. The 
broader dispersion suggests that the Si IV may be 
in slightly lower density gas than the Si II, since the 
pressure is presumably the same for both. It is not 
clear however that the observed difference in velocity 
width between Si IV and Si II is real. One problem 
is that the continuum around Si IV appears less well 
defined than for the Si II and Si III features. As else- 
where in this paper, we assume a stellar continuum 
which is linear in logF-logA, but in fact there is a 
hint of curvature, a large scale depression in the con- 
tinuum around Si IV 1394, 1403 A (see WCHFLS96, 
Figure 1). If we have systematically misjudged the 
continuum, then it is possible that we have underes- 
timated the uncertainties in the parameters of Si IV 
given Table |^, perhaps by as much as a factor of 2. 
The Gaussian profile shown in Figure ^ is constrained 
to have the same center and 1240 kms^^ dispersion 
as the Si II 1260 A feature. Visually, at least, the fit 
appears satisfactory. 

2.7. Total Si mass 



In subsection 2.5 we showed that the observed col- 



umn density of shocked Si II is close to the theoretical 
steady state value. Is the same also true for Si III and 
Si IV? The answer is no. In steady state, the pre- 
dicted column densities of shocked ionic species are 
inversely proportional to their respective coUisional 



ionization rates {av) , modified by an appropriate elec- 
tron to ion ratio (cf. equation ||l3]). Adopting rates 
6.1xl0-8cm3s"\ 2.2xl0~^cm3", 1.1 x lO^^cm^s"! 
for Si II, Si III, Si IV respectively (Lennon et al. 1988), 
and assuming a nominal i electrons per ion for Si^% 
yields relative column densities in steady state 



Nsiu ■ Nsiiu : A^siiv 



1 : 1.4 : 1. 



(15) 



By comparison the observed shocked column densities 
are A^siii : ^Siiii : A^Siiv = 1 : 0.43 : 0.41, according 
to Table ^ Evidently the observed abundances of 
shocked Si III and Si IV relative to Si II are several 
times less than predicted in steady state. 

As discussed in subsection 2.5, there are several 



ways to reduce the column density below the steady 
state value, of which the most obvious is to truncate 
the column density, as is strongly suggested by the 
fact that the ionization timescales of Si III and Si IV 
are becoming long compared to the age of the rem- 
nant. In fact these ionization timescales are in pre- 
cisely the same ratio (cf. equation [Q) as the steady 
state column densities (nsh 



(16) 



isiii ■ ^Siiii : ^Siiv = 1 : 1.4 : l.i 



and it has already been seen that the ionization 
timescale tsni of Si II is comparable to the age of 
the remnant, equation (^). 

If it is true that it is the long ionization times 
which cause the column densities of Si III and Si IV 
to be lower than steady state, then this suggests that 
there may be little Si in higher ionization states in 
the shocked gas on the far side of SN1006. Thus it 
appears plausible that we are observing in UV ab- 
sorption essentially all the Si there is on the far side 
of SN1006 along the line of sight to the background 
SM star. To avoid confusion, it should be understood 
that this statement refers specifically to Si on the far 
side along this particular line of sight. Higher ion- 
ization states of Si are indicated by the observed Si 
x-ray line emission (Koyama et al. 1995), which could 
arise from denser shocked gas in other parts of the 
remnant. 

The mass of Si can be inferred from the observed 
column densities of Si II, Si III, and Si IV, if it is 
assumed that silicon was ejected spherically symmet- 
rically by the supernova explosion. We will argue in 
subsection 2.8 that the absence of blueshifted absorb- 



ing Si is not inconsistent with spherical symmetry. 
Taking the shocked and unshocked masses of Si II 
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Table 2 

Column densities of shocked and unshocked Si III and Si IV, relative to best fit column densities 

OF Si II 





Si II 


Si III 


Si IV 


Shocked 


1±0.03 


0.43 ±0.02 


0.41 ±0.02 


Unshocked 


1 ±0.13 


0.065 ±0.035 


0.02 ±0.07 



Note. — Absolute column densities of Si II are 
given in Table |l|. Column densities of Si IV rela- 
tive to Si II are for fits in which the dispersion of 
Si IV is constrained to be that of Si II 1260 A, namely 
1240 kms~^. The column densities of Si IV become 
0.46 ± 0.02 (shocked) and -0.06 ± 0.06 (unshocked) if 
instead the dispersion of the Si IV is taken to have the 
best fit value 1700 ± 100 kms"^ 



from Table 
from Table 



and the ratios of Si II, Si III, and Si IV 
yields a total inferred Si mass of 



Msi = 0.127(1 + 0.43 + 0.41)Mq 
= 0.25 ± 0.01 Mq . 



0.017(1 



O.O65)M0 
(17) 



This is comparable to the 0.16 Mq of Si in model W7 
of Nomoto et al. (1984). It is also consistent with 
the 0.20 Mq of Si inferred from the strength of the 
SiK line observed with ASCA (Koyama et al. 1995). 
Koyama et al. do not quote a mass, but they do state 
that the measured surface brightness of the Si K line 
is 5 times higher than that in the modelQ of Hamilton, 
Sarazin & Szymkowiak (1986; see also Hamilton et al. 
1985), which had 0.04 Mq of Si. 



2.8. No blueshifted Si 

homogeneous ISM? 



evidence for an in- 



There is no evidence for blueshifted Si absorp- 
tion in the UV spectrum. At best, there is a pos- 
sible hint of a broad shallow depression around ~ 



^ Hamilton et al. took the strength of the Si K line to be (the 
upper limit to) that measured from the Einstein Solid State 
Spectrometer (Becker et al. 1980), so there is a discrepancy 
between the ASCA and SSS data. However, Hamilton et al. 
also noted in their Figure 6 a marked discrepancy between the 
SSS and HEAO-1 data of Galas, Venkatesan & Garmire (1982), 
so it is reasonable to suspect an error in the normalization of 
the SSS data. 



-5000kms-\ from 1365 A to 1390 A, on the blue 
side of the Si IV 1394, 1403 A line (see WCHFLS96, 
Figure 1). The possibility that some high velocity 
blueshifted Si II 1260 A is hidden in the red wing of 
Si III 1206 A is excluded by the absence of correspond- 
ing blueshifted Si II 1527 A. 

There are two possible reasons for the asymmetry 
in the observed Si absorption. One is that there was 
an intrinsic asymmetry in the supernova explosion. 
According to Garcia-Senz & Woosley (1995), the nu- 
clear runaway that culminates in the explosion of a 
nearly Chandrasekhar mass white dwarf begins as sta- 
ble convective carbon burning, and ignition is likely to 
occur off-center at one or more points over a volume of 
the order of a pressure scale height. The subsequent 
propagation of the convectively driven burning front 
is Rayleigh- Taylor unstable (Livne 1993; Khokhlov 
1995; Niemeyer & Hillebrandt 1995), although Ar- 
nett & Livne (1994) find that in delayed detonation 
models the second, detonation, phase of the explosion 
tends to restore spherical symmetry. 

Thus asymmetry in the explosion is possible, per- 
haps even likely, as an explanation of the asymmetry 
in the Si absorption, especially since the absorption 
samples only a narrow line of sight through SN1006 to 
the background SM star. However, we do not pursue 
this possibility further, in part because in abandoning 
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spherical symmetry we lose any predictive power, and 
in part because there is another explanation which is 
more attractive because it resolves some other obser- 
vational problems. 

The other possible cause of the asymmetry in the 
Si absorption is that the ISM around SN1006 is inho- 
mogeneous, with the ISM on the far side of SN1006 
having a significantly lower density than the near side. 
According to this hypothesis, the low density on the 
far side is such that the reverse shock on the far side 
has reached inward only to a free expansion radius of 
7070 km s~^, whereas the higher density on the near 
side is such that the reverse shock on the near side 
has passed entirely through the Si layer, and Si has 
been collisionally ionized to stages higher than Si IV, 
making it unobservable in UV absorption. 

A serious objection to the reverse shock on the 
near side being farther in than the reverse shock 
on the far side is the observation by WCFIIS93 of 
blueshifted Fe II absorption certainly to velocities 
— 7000kms~^, perhaps to velocities — 9000kms~^. 
We will review the observational evidence for such 
high velocity blueshifted Fe II in Section ^ below, 
where we will argue that the evidence is not com- 
pelling. 

An inhomogeneous ISM around SN1006 is indi- 
cated by other observations. Wu et al. (1983) and 
subsequent authors have remarked on the difhculty 
of fitting the observed high velocity (~ 7000 kms~^) 
Si within the confines of the observed interstellar 
blast wave, if spherical symmetry is assumed. At a 
distance of 1.8 ± 0.3 kpc (Laming et al. 1996), the 
remnant's observed 15' radius (Reynolds & Gilmore 
1986, 1993; see also Wilhngale et al. 1996) at 980 
years oldQ corresponds to a free expansion radius of 
7700 ± 1300 km s-^ The difficulty is resolved if the 
remnant of SN1006 bulges out on the far side, because 
of the lower density there. 

A second piece of evidence suggesting inhomogene- 
ity is the high 5050 km s~^ velocity of the shocked 
ejecta behind the reverse shock on the far side of 
SN1006 inferred from the present observations, com- 
pared with the 1800-2400 kms^^ (the lower velocity 
corresponds to no coUisionless electron heating, which 
is the preferred case) velocity of shocked gas behind 

^ Note there is a 23 year light travel time across one radius of 
the remnant, so really we are seeing the far side at an age 23 
years younger, and the near side 23 years older, than the mean 
age. 



the interstellar shock inferred from H a and UV emis- 
sion line widths along the NW sector of the remnant 
(Kirshner et al. 1987; Long et al. 1988; Smith et al. 
1991; Raymond et al. 1995). These two velocities, 
5050 km s~^ versus 1800-2400 kms~^, appear incom- 
patible, especially as the velocity of shocked gas is ex- 
pected in realistic hydrodynamic models to increase 
outwards from the reverse shock to the interstellar 
shock (see for example HF88, Fig. 2). The incompat- 
ibility is resolved if the ISM around SN1006 is inho- 
mogeneous, with the density on the far side of SN1006 
being substantially lower than the density at the NW 
edge. 

A final piece of evidence supporting inhomogene- 
ity comes from the observation of SiK line emission 
in x-rays (Koyama et al. 1995). This emission is likely 
to be from ejecta, since the inferred abundance of Si 
is stated to be an order of magnitude higher than 
that of O, Ne, or Fe. To ionize Si to high ioniza- 
tion states in the age of SN1006, and to produce SiK 
line emission at the observed luminosity, requires den- 
sities nsi ^ lO^^cm"'^ (i.e. rig ^ lO^^cm"^ since 
the Si is highly ionized), substantially higher than 
the postshock density of ngin — 2 x lO^'' cm~'^ de- 
duced here from the redshifted Si II 1260 A absorp- 
tion profile. It is difficult to see how the required high 
Si density could be achieved if the reverse shock ev- 
erywhere around the remnant is at a radius as high 
as 7070 km s^^, whereas higher densities would oc- 
cur naturally if the reverse shock around most of the 
remnant were farther in, since then ejecta would have 
been shocked at earlier times when its density (oc 
in free expansion) was higher. 

All these arguments point to the notion that the 
ISM density on the far side of SN1006 is anomalously 
low compared to the density around the rest of the 
remnant. 

2.9. Ionization of Si on the near side 

In Section |^, we will argue that the Fe II absorp- 
tion profiles suggest that the reverse shock on the near 
side of SN1006 may be at a free expansion radius of 
4200 kms^^ (as with the Hubble expansion of the Uni- 
verse, it is often convenient to think in a comoving 
frame expanding freely with the ejecta, so that a free 
expansion velocity v = r/t can be thought of as a ra- 
dius) . Here we estimate whether Si could have ionized 
beyond Si IV, as required by the absence of observed 
blueshifted Si absorption, if the reverse shock on the 
near side is indeed at 4200 km s~^. 
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In our first estimate, we find that blueshifted Si IV 
should be observable, with a column density ~ 40% 
that of the observed redshifted Si IV absorption. 
However, the column density is somewhat sensitive 
to the assumptions made, and it is not hard to bring 
the column density down below observable levels. 

The line profile of any blueshifted absorbing Si 
would have a width comparable to that of the ob- 
served broad redshifted Si absorption, but the cen- 
troid would be shifted to lower velocities, to ^ —2000 kms~^, 
if it assumed that the reverse shock velocity on the 
near side is comparable to that, Vg = 2860 km s""'^, 
observed on the far side. To avoid being observed at 
say 3(7, the column density of blueshifted Si IV should 
be less than 0.6 x 10^"' cm^^. 

To estimate the ionization of Si, it is necessary to 
adopt a hydrodynamic model. Now one interesting 
aspect of hydrodynamical models of defiagrated white 
dwarfs expanding into a uniform ambient medium is 
that the reverse shock velocity remains almost con- 
stant in time (this conclusion is based on hydrody- 
namic simulations carried out by HF88). In model 
W7 of Nomoto et al. (1984), the reverse shock veloc- 
ity varies (non-monotonically) between 3300 km s~^ 
and 5200 km s^^ as it propagates inward from a free 
expansion radius of 13000 kms~^ to 700kms~^, af- 
ter which the shock velocity accelerates. Similarly in 
model CDTG7 (Woosley 1987, private communica- 
tion), which is similar to model CDTG5 of Woosley 
& Weaver (1987), the reverse shock velocity varies be- 
tween 3200 kms~^ and 4100 kms^^ as it propagates 
inward from a free expansion radius of 10000 kms~^ 
to 1500 km s~^. These numbers do not depend on the 
density of the ambient medium, although they do de- 
pend on the ambient density being uniform. If the 
reverse shock velocity Vs remains constant in time, 
then the radius r of the reverse shock evolves with 
time t according to dr/dt = r/t — Vs, from which it 
follows that the free expansion radius r/toi the re- 
verse shock varies with time as 



— — Vs hileft{-^right) 



(18) 



where is the age at which the reverse shock even- 
tually hits the center of the remnant. 

The assumption that the reverse shock ti^ is con- 
stant in time may not be correct for SN1006, but it 
provides a convenient simplification to estimate the 
ionization of Si on the near side. Let us first esti- 
mate the ionization state of Si which was originally 



at a free expansion radius of 7070 km s~^, the current 
location of the reverse shock on the far side. If the 
reverse shock on the near side is currently at a free 
expansion radius of r/t = 4200 kms^^, and if it has 



moved at a constant = 2860 km s 



the measured 



value on the far side (recall that the reverse shock 
velocity is independent of the ambient density, for 
uniform ambient density), then it would have passed 
through a free expansion radius of 7070 km s~^ when 
the age of SN1006 was t^/t = exp[(4200kms"i - 
7070kms"i)/2860kms"i] = 0.37 times its present 
age t, according to equation (|l8|). The postshock 
density of Si ions at that time would have been 
(t/ts)^ = 20 times higher than the presently observed 
postshock density of nsiii = 2.2 x 10"** cm^"^. The ion 
density in the parcel of gas shocked at that time has 
been decreasing because of adiabatic expansion. The 
rate of decrease of density by adiabatic expansion can 
be inferred from the observed global rate of expansion 
of the remnant (Moffet, Goss & Reynolds 1993) 



i?oci" , a = 0.48 ±0.13 



(19) 



which for an assumed uniform ambient density would 
imply that the pressure is decreasing as P cx (R/t)'^ oc 
^2a-2^ hence that the density in Lagrangian gas ele- 
ments is varying as n oc p3/5 ^ ^6(a-i)/5^ ^pj^^ 
rent ionization time r = ndt oi the parcel of gas 
originally at free expansion radius 7070 km s~^ which 
was shocked at an age tg/t — 0.37 is then 



"Sill* 



(6a- l)/5 \t 
= 6.0 X lO'^cm^^s 



(6q-1)/5 



(20) 



where nsiiii = 6.6 x 10® cm^^ s is the present post- 
shock density of Si II ions times age at the radius 
7070 km s~^. If the Si is assumed unmixed with other 
elements and initially singly ionized, then the ioniza- 
tion time d^) yields ion fractions 5% Si III, 34% Si IV, 
and the remaining 61% in higher stages. 

Since this ionization state is close to (just past) 
the peak in the Si IV fraction, it should be a good 
approximation to estimate the column density of 
Si IV by expanding locally about the conditions at 
7070 kms~^. The expected column density is the 
steady state column density, multiplied by a geomet- 
ric factor, and further multiplied by a 'density pro- 
file' factor [1 -I- / {nsts)]^^ , as shown in the Ap- 
pendix, equation (A7). The steady column density 
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steady 

SilV — 



nsiiv'fs/(4ne(CTw)siiv) = 22 x 10^ 



is N^^^^ 

assuming a nominal n-e/nsiw — 3. The geomet- 
ric factor is (7070^^/4200^)2 = 0.38, which is the 
squared ratio of the radius of the gas at the time 
ts it was shocked to its radius at the present time 
t. For the density profile factor, equation (^) gives 
r/(n,t,) = [it/U)('^'--^y^ - l]/[(6a - l)/5] = 1.22, 
while the logarithmic slope 7 = —d\nn/d\nts\t of 
the shocked density profile, equation (A3), is 7 — 
3+t;^ainn™"h/5(r/t)+6(a-l)/5 = 3.7, the 3 coming 
from free expansion, the Vsd\nn^'^^^/d{r/t) = 1.37 
from the observed unshocked Si II density profile (see 
eq. (2^) at 7070 kms~^ along with equation (18), and 
the 6(0; — — —0.62 from the reduction in density 
caused by adiabatic expansion. The resulting den- 
sity profile factor is [1 + jT/{nsts)]^^ — 0.18. The 
expected column density of blueshifted Si IV is then 
iVgiiv = 1.5 X 10"^'* cm~2, which is 40% of the observed 
column density of redshifted Si IV, and 2.5 times the 
minimum (3(t) observable column density. 

Thus under a reasonable set of simplifying assump- 
tions, there should be an observable column density of 
blueshifted Si IV, contrary to observation. However, 
the expected column density is sensitive to the as- 
sumed parameters. For example, reducing the shock 
velocity on the near side by 20% to 2300 kms""'^ re- 
duces the column density by a factor 2.5 to the ob- 
servable limit A^siiv — 0.6 x 10^** cm^^. Whether the 
shock velocity on the near side is less or more than 
that on the far side depends on the unshocked density 
profile of ejecta (generally, a shorter exponential scale 
length of unshocked density with velocity yields lower 
shock velocities). The expected column density is also 
sensitive to the shocked density profile, as might be 
guessed from the fact that the density profile factor of 
0.18 estimated above differs substantially from unity. 

Alternatively, the column density of Si IV could 
be reduced below the observable limit by mixing the 
Si with a comparable mass of other elements, such 
as iron, since this would increase the number of elec- 
trons per silicon ion, causing more rapid ionization. 
The possibility that there is iron mixed with Si at ve- 
locities ^ 7070 kms~^ gains some support from the 
observed density profile of Fe II, shown in Figure ^ 
below. Note this does not conflict with the argument 
in subsection 2.5 that most of the shocked Si (which 
was originally at higher free expansion velocities) is 
probably fairly pure, with little admixture of other 
elements such as iron. 



3. Iron 

We have argued above that an attractive explana- 
tion for the presence of redshifted Si absorption and 
absence of blueshifted absorption is that the ISM on 
the near side of SN1006 is much denser than that 
on the far side, so that a reverse shock has already 
passed all the way through the Si layer on the near 
side, ionizing it to high ionization stages, whereas the 
reverse shock is still moving through the Si layer on 
the far side. This picture appears to conflict with 
our previously reported result (WCFHS93), accord- 
ing to which blueshifted Fe II is present to velocities 

— 8000kms^^. In this Section we reexamine the 
Fe II absorption lines to see how robust is this result. 

In the average broad Fe II profile shown in Fig- 
ure 3 of WCFHS93, redshifted Fe II seems to extend 
up to about 7000 km s~^, but not much farther. This 
is consistent with the argument of the present pa- 
per, which is that the reverse shock on the far side 
of SN1006 lies at 7070 kms'^ The problem lies on 
the blueshifted side of the Fe II profile, which ap- 
pears to extend clearly to — 7000kms~^, possibly to 
-9000kms-i. 

Figure 2 of WCFHS93 shows separately the two 
broad Fe II 2383, 2344, 2374 A and Fe II 2600, 2587 A 
features. The Fe II 2600, 2587 A feature appears to 
have a fairly abrupt blue edge, although there is per- 
haps a tail to higher velocities depending on where the 
continuum is placed. The blue edge is at a velocity 
of — 4200 km s^^ with respect to the stronger 2600 A 
component of the doublet, and the same edge appears 
at this velocity in the average Fe II profile shown in 
WCFHS93's Figure 3. We will argue that this edge 
plausibly represents the position of the reverse shock 
on the near side of SN1006. 

In contrast to Fe II 2600 A, the deconvolved pro- 
file of the Fe II 2383 A feature, plotted in the bottom 
curve of WCFHS93's Figure 2, shows blueshifted ab- 
sorption clearly extending to 7000 km s^^. We note 
that the second strongest component of the triplet, 
2344 A, with 1 /3 the oscillator strength of the princi- 
pal 2383 A component, lies at — 4900kms~^ blueward 
of the principal line, and uncertainty involved in re- 
moving the secondary component in the deconvolu- 
tion procedure could tend to obscure any sharp blue 
edge at — 4200kms"^ on the principal component. 



14 




0x10-14 n 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 n 

2000 2100 2200 2300 2400 2500 2600 2700 2800 2900 3000 

Wavelength (A) 

Fig. 4. — G190H and G270H spectra, with the cahbratcd spectra at bottom, and the dereddened spectra at 
top. Also shown are the continua Unear in logF-logA (sohd hnes) adopted in the present paper for each of 
the Fe II 2600 A and Fe II 2383 A features, and a continuum quadratic in F-X (dashed Une) similar (but not 
identical, because of the slightly different reddening) to that adopted by WCFHS93. The linear continua are 
logF = log(1.89 X 10-i4ergs-icm-2 A"^)-2.71og(A/2383A) for the Fell 2383 A feature, andlogF = log(1.38 x 
10-"ergs-i cm-2 A"^) - 3.65 log(A/2600 A) for the Fe II 2600 A feature. 



3.1. Fe II analysis 

In this subsection we present details of a reanal- 
ysis of the Fe II absorption Hnes in the HST G190H 
and G270H spectra originally analyzed by WCFHS93. 
The observations are described by WCFHS93, and 
here we describe the differences between the analysis 
here and that of WCFHS93. In carrying out the re- 
analysis, we paid particular attention to procedures 
which might affect the blue wing of the Fe II 2383 A 
feature. 

The G190H and G270H spectra overlap over the 

wavelength range 2222-2330 A, and we merged the 
two spectra in this region using inverse variance 
weighting, whereas WCFHS93 chose to abut the spec- 
tra at 2277 A. In merging the spectra we interpolated 
the G270H data to the same bin size as the G190 spec- 
trum, which has higher resolution (2 A versus 2.8 A), 
and higher signal to noise ratio than the G270H spec- 
trum in the overlap region. According to the FOS 



handbook, there is contamination at the level of a few 
percent in the G190H spectrum above 2300 A from 
second order, but, given the absence of strong fea- 
tures over 1150-1165 A, we accept this contamination 
in the interests of obtaining higher signal to noise ra- 
tio. The merged spectrum is noticeably less choppy 
than the G270H spectrum alone in the overlap region. 

The 2200 A extinction bump is close to the blue 
wmg of the Fe II 2383 A feature, so we re-examined 
the reddening correction. In practice, the changes 
made here had little effect on the profile of the 
Fe II 2383 A feature. We dereddened the G190H and 
G270H spectra using the extinction curve of Cardelli, 
Clayton & Mathis (1989), adopting Eb-v = 0.119, 
which is the best fitting value determined by Blair et 
al. (1996) from HUT data, and R = Ay /Eb-v = 3.0. 
The value of Eb-v if' slightly higher than the value 
Eb-v = 0.1 ± 0.02 measured by WCFHS93 us- 
ing the extinction curve of Savage & Mathis (1979). 
WCFHS93 comment that their dereddening leaves a 
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bump from 1900 A to 2100 A and a shallow trough 
from 2100 to 2300 A. We find the same difficulty here: 
the slightly higher value of Eb-v adopted here does 
slightly better at removing the 2200 A depression, but 
at the expense of producing a bigger bump at 2000 A. 
The choice of R makes little difference, but lower val- 
ues help to reduce the bump marginally. The value 
R = 3.0 adopted here is slightly below the 'standard' 
value R = 3.1. 

WCFHS93 fitted the continuum flux to a quadratic 
function of wavelength. The simple quadratic form 
does impressively well in fitting the dereddened spec- 
trum over the full range 1600 A to 3300 A (cf. Figure | 
and WCFHS93 Figure 2). However, the quadratic 
form does not fit perfectly, and there remains a resid- 
ual discrepancy which is not well fitted by a low order 
polynomial, and which may possibly result from im- 
perfect modeling of the extinction, especially around 
the 2200 A bump. The imperfection mainly affects 
the Fe II 2383 A feature: the quadratic contimmm 
appears too steep compared to the 'true' continuum 
around this feature. Here we resolve the difficulty 
by the expedient of fitting the dereddened continua 
around each of the two broad Fe II features sepa- 
rately, to two different linear functions in logF-logA. 
The adopted continua are shown in Figure ^ 

An important difference between the present anal- 
ysis and that of WCFHS93 is in the treatment of 
narrow lines. WCFHS93's procedure was to iden- 
tify all lines with an equivalent width, defined rela- 
tive to a local continuum, greater than 3 times the 
expected error. WCFHS93 then subtracted the best 
fitting Gaussian for each such line, treating the po- 
sition, width, and strength of each line as free pa- 
rameters. Here we adopt a different policy, requiring 
that the positions, widths, and strengths of narrow 
lines conform to prior expectation. That is, for each 
identified narrow line we subtract a Gaussian profile 
in which the wavelength is set equal to the expected 
wavelength (modulo an overall -1-36 km s"^ shift for all 
lines), the width is set equal to the instrumental res- 
olution (2.8 A FWHM for G270H), and the strengths 
are required to be mutually consistent with other nar- 
row lines of the same ion. 

The relevant narrow lines are those in and near the 
broad Fe II features. In the Fe II 2383, 2344, 2374 A 
feature, besides the narrow (presumed interstellar) 
components of the Fe II lines themselves, we identify 
the narrow fine at 2298 A as steUar C III 2297.578 A 
(Bruhweiler, Kondo & McCluskey 1981). WCFHS93 



subtracted an unidentified narrow line at 2316 A, but 
the G190H spectrum does not confirm the reality of 
this line in the G270H spectrum, and here we leave 
it unsubtracted. In the Fe II 2600, 2587 A feature, 
besides the narrow Fe II lines themselves, we iden- 
tify narrow fines of Mn II 2577, 2594, 2606 A, as did 
WCFHS93. 

The mean velocity shift of the three most promi- 
nent narrow Fe II lines, those at 2383 A, 2600 A, and 
2587 A, is +36 ± 24kms^^, and we adopt this ve- 
locity shift for all the narrow Fe II and Mn II lines. 
We aUow the stellar C III 2298 A line its own best fit 
velocity shift of +19kms~^, since there is no reason 
to assume that the stellar and interstellar velocities 
coincide exactly. 

The observed equivalent widths of the Mn II lines 
are approximately proportional to their oscillator strengths 
times wavelengths, which suggests the lines are unsat- 
urated, so we fix the ratios of the fitted Mn II lines at 
their unsaturated values. 

Of the five narrow Fe II lines, the two principal 
lines Fe II 2383 A and Fe II 2600 A, and also the 
line with the fourth largest oscillator strength, Fe II 
2587 A, have approximately equal observed equiv- 
alent widths (in velocity units) relative to a local 
continuum, Vl^2383, W^2600, W^2587 = 85 ± 6kms~\ 
64 ± 6kms~^, 62 ± 6kms~^ respectively (at fixed 
centroid and dispersion), which suggests the lines are 
saturated. The fifth and weakest line, Fe II 2374 A, 
has an observed equivalent width about half that of 
the strong lines, which is consistent with the weak 
line being marginally optically thin and the strong 
lines again being saturated. For these four lines we 
allow the strength of the fitted line to take its best fit 
value, since they are mutually consistent within the 
uncertainties. 

The line with the third largest oscillator strength, 
Fe II 2344 A, appears anomalous, since the observed 
line has an equivalent width less than 1/4 that of the 
strong lines, or 1/2 that of the intrinsically weaker 
Fe II 2374 A line. In the fit, we force the equiva- 
lent width of the anomalous Fe II 2344 A narrow line 
to be the saturated value measured from the Fe II 
2600 A and Fe II 2587 A fines, multiplied by 0.8 to 
allow for a 2a uncertainty in this value. The fit gives 
the impression that the anomalous Fe II 2344 A line 
is oversubtracted, but the effect is to bring the profile 
of the deconvolved broad Fe II 2383 A line into closer 
agreement with that of Fe II 2600 A. 
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Fig. 5.— G270H spectrum showing at left the Fe II 2382.765, 2344.214, 2374.4612 A (/ = 0.3006, 0.1097, 0.02818; 
we ignore an unobserved fourth component Fe II 2367.5905 A with / = 1.6 x 10^^) feature, and at right the Fe II 
2600.1729, 2586.6500A (/ = 0.2239, 0.06457) feature (Morton 1991). Below 2330 A, the spectrum is an inverse- 
variance-weighted merger of G19011 and G270H spectra. The lower curve shows the dereddened spectrum, scaled 
to a continuum, with narrow interstellar Fe II, stellar C III 2297.578 A, and interstellar Mn II 2576.877, 2594.499, 
2606.462 A (/ = 0.3508, 0.2710, 0.1927) lines subtracted as indicated by dashed lines. The upper curve (offset 
by log 1.15) shows the deconvolved spectra, after removal of the weaker components of the broad Fe II lines. The 
velocity scales shown on the upper axis are with respect to the rest frames of the principal component of each of 
the features, the Fe II 2383 A line on the left, and the Fe II 2600 A line on the right. The adopted continua are as 
shown in Figure ^. 



We deconvolved the broad Fe II 2383, 2344, 2374 A 
and Fe II 2600, 2587 A features by subtracting the 
contributions from the weaker components, using the 
following analytic procedure. In a two component 
line, the observed optical depth t{v) at velocity v with 
respect to the principal component is a sum of the line 
profile </>(v) of the principal component and the line 
profile e4>{v + Av) of the secondary component, where 
Av is the velocity shift of the secondary relative to the 
principal component, and e = f2^2/ (fi^i) < 1 is the 
ratio of oscillator strengths times wavelengths: 



t{v) = 4>{v) + e0(i' + Aw) 



(21) 



Equation ( |2l| ) inverts to 

(j)[v) = t{v) - tT{v + Av) + t^T{v + 2Av) H (22) 

which can conveniently be solved iteratively by 



(j)i(v) 



'{v) — er(t) + Av) 



,2" 



Mv + 2"Av) . (23) 



The iterative procedure converges rapidly to the solu- 
tion, (pn 4>, as n increases; we stop at (j)'^. To avoid 
irrelevant parts of the spectrum outside the line pro- 
file from propagating through the solution, we set the 
optical depth to zero, riy) = in equation (p3|), at 
velocities v > 10000 km s^^. The above procedure 
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works for a two component line such as Fe II 2600, 
2587 A, and a slightly more complicated generaliza- 
tion works for a three component line such as Fe II 
2383, 2344, 2374 A. 

3.2. Fe II line profiles 

Figure || shows the results of our reanalysis. The 
upper curves in the Figure show the deconvolved Fe II 
2383 A and Fe II 2600 A line profiles, and these decon- 
volved profiles agree well with each other. The decon- 
volved Fe II 2600 A profile here also agrees well with 
that of WCFHS93. However, the revised Fe II 2383 A 
profile no longer shows compelling evidence for high 
velocity blueshifted absorption beyond —4500 kms"^, 
although the presence of some absorption is not ex- 
cluded. 

What causes the difference between the Fe II 
2383 A line profile shown in Figure || versus that of 
WCFHS93? One factor is that we adopt different con- 
tinua, as illustrated in Figure H. WCFHS93's single 
quadratic fit to the continuum over the entire spec- 
trum is certainly more elegant than the two separate 
linear fits to each of the two broad Fe II features which 
we adopt here. The advantage of the fit here is that 
it removes the apparent tilt in the Fe II 2383 A line 
profile left by the quadratic fit, evident in Figures 2 
and 3 of WCFHS93. 

However, the major difference between the two 
analyses is the subtraction here of the narrow Fe II 
2344 A interstellar line with a strength 0.8 times the 
saturated line strength observed in the Fe II 2600 A 
and 2587 A narrow lines. By comparison, WCFHS93 
subtracted the narrow Fe II 2344 A line using the ob- 
served strength of the line, which is anomalously weak 
compared to the other four narrow Fe II lines. 

Figure ^ shows the Fe II density inferred from the 
mean of the two deconvolved Fe II features. The Fe II 
column density inferred from the mean profile, inte- 
grated from — 4500kms^^ to +7100 kms^^, is 



TVfoII = 10.8 ± 0.9 X 10^"* cm 



(24) 



Most of the uncertainty, based here on the scatter 
between the two deconvolved profiles, comes from the 
blue side of the profile: the column density on the 
blue side from — 4500 km s^^ to Okms~^ is A^foII = 
5.2 ± 0.8 X lO^^cm"^, while the column density on 
the red side from Okms^^ to +7100 kms~^ is TVpcii = 
5.6 + 0.2 X 10^''cm-2. 

In estimating the mass of Fe II from the density 



profile in Figure ^ we take into account the small 
correction which results from the fact that the SM 
star is offset by 2'. 45 ± 0'.25 southward from the pro- 
jected center of the 15' radius remnant (Schweizer & 
Middleditch 1980). The offset corresponds to a free 
expansion velocity of v± ~ 1300 ± 250kms~^ at the 
1.8 ± 0.3 kpc distance of the remnant. If spherical 
symmetry is assumed, then the mass is an integral 
over the density p{v) at line-of-sight velocity v. 

(■-"max 

M = M(<vi_)+t^ p{v){v^ + vl)^^^ vdv (25) 



where M(<v±) is the mass inside the free-expansion 
velocity v±. At a constant central density of ppoii = 
0.005 X 10~^^ gmcm"'^, the mass inside i;^ = 1300 kms~^ 
would be MfcIi(< v±) = 0.0007 M©, and the ac- 
tual Fe II mass is probably slightly higher, given 
that the density is increasing mildly inward. The 
masses given below, equation (^6|), include a fixed 
MFeii(< vj_) = 0.001 M0. The factor {v^ + )i/2 



rather than v in equation (25) increases the Fe II 
masses by a further 0.002 Mq, so the masses quoted 
in equation ( |2^ ) are altogether 0.003 Mq larger than 
they would be if no adjustment for the offset of the 
SM star were applied. 

The total mass of Fe II inferred from the cleaner, 
red side of the profile, assuming spherical symmetry, 
is then 



Ah 



0.0156 + 0.0009 Mq (w < 4200kms-i) 
Fell = <v 0.0195 + 0.0013 Mq (w < 5000kms-i) 
0.029 + 0.004Mq (w < 7070kms-i) . 

(26) 

The uncertainties here are based on the scatter be- 
tween the two deconvolved Fe II profiles, and do not 
included systematic uncertainties arising from place- 
ment of the continuum, which mostly affects the 
outer, high velocity parts of the profile. WCFHS93 
obtained MfeIi = 0.014 Mq from the red side of the 
mean Fe II profile, which is lower than the A/fcIi = 
0.029 Mq obtained here mainly because of the differ- 
ent placement of the continuum (see Figure |4|) , and 
to a small degree because of the adjustment applied 
here for the offset of the SM star. 

Figure ^ also shows for comparison the density of 
unshocked Si. Below 7070 kms~^, the Si density pro- 
file is just the unshocked profile inferred from the Si II 
1260 A absorption. Figure |l|. Above 7070kms"\ the 
Si density is a plausible extrapolation which is consis- 
tent with observational constraints: it is a quadratic 
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Table 3 

Parameters measured from Fe II profile 



Parameter 


Value 


Position of reverse shock on near side 


-4200± lOOkms-i 


Column density of Fe II 


10.8 ±0.9 X lOi^'cm-^ 


Mass of Fe II up to 7070 kms^^ 


0.029 ± 0.004 Mq 



Note. — Fe II mass is from red side of profile, and assumes 
spherical symmetry. 



function of the free expansion velocity v 



0.00413 X lO^^'^gmcm"^ 



{v - 5600kms~^)(12000kms~^ - v\ 



(3200 kms-i)2 



(27) 



which approximately reproduces the observed profile 
of unshocked Si II, and which contains, on the as- 
sumption of spherical symmetry, a total Si mass of 
0.25 Mq, in accordance with equation (p^). 

3.3. Reverse shock on the near side 

If the reanalysis of the Fe II lines here is accepted, 
then it is natural to interpret the sharp blue edge on 
the broad Fe II lines at —4200 km s~^ as representing 
the free expansion radius of the reverse shock on the 
near side of SN1006. This identification is not as con- 
vincing as the identification of the sharp red edge on 
the Si II 1260 A feature as representing the radius of 
the reverse shock on the far side at 7070 km s~^. 

Figure illustrates schematically the inferred struc- 
ture of the remnant of SN1006. The picture is in- 
tended to be approximately to scale, and in it the 
diameter of SN1006 along the line of sight is roughly 
20% larger than the diameter transverse to the line of 
sight. By comparison, the diameter of the observed 
radio and x-ray remnant varies by 10%, from a mini- 
mum of 30' to a maximum of 33' (Reynolds & Gilmore 
1986, 1993) or 34' (Wilhngale et al. 1996). 



As already discussed in subsections and p.9| , 
if the position of the reverse shock on the near side 
at 4200 kms^^ is typical of the rest of the remnant, 
while the 7070 km s^^ position of the reverse shock 



on the far side is anomalously high because the in- 
terstellar density on the far side is low, then it would 
resolve several observational puzzles. In summary, 
these observational puzzles are: (1) how to fit gas ex- 
panding at 7000 kms""'^ within the confines of the 
interstellar shock (answer: the remnant bulges out 
on the far side because of the low density); (2) how 
the 5050 kms~^ velocity of shocked Si on the far side 
could be so much higher than the 1800 kms^^ ve- 
locity (assuming no coUisionless electron heating) of 
gas behind the interstellar shock along the NW fila- 
ment (answer: velocities on the far side are anoma- 
lously high because the interstellar density there is 
anomalously low); (3) how to achieve Si densities 
"Si ^ 10~^ cm~^ necessary to produce the observed Si 
x-ray emission, compared to the postshock Si density 
of 2.2 X 10~'*cm~^ measured from the Si II absorp- 
tion on the far side (answer: gas shocked at earlier 
times is denser because density decreases as p cx t'^ 
in free expansion); and (4) why there is strong red- 
shifted Si absorption but no blueshifted absorption 
(answer: Si on the near side has been shocked and 
coUisionally ionized above Si IV) . As regards the sec- 
ond of these problems, if the reverse shock on the 
near side is indeed at 4200 km s~^, then the veloc- 
ity of reverse-shocked gas on the near side would be 
of order 2000 km s~^, much more in keeping with the 
1800 kms^^ velocity of shocked gas in the NW fila- 
ment. 

3.4. Contribution of shocked Fe to absorption 



In subsections 2.1 and 2.6 we concluded that most 



of the observed absorption by Si ions is from shocked 
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Si. Does shocked Fe also contribute to the observed 
broad Fe II absorption profiles? The answer, on both 
observational and theoretical grounds, is probably not 
much. 

On the red side of the Fe II profile, shocked 
Fe II would have a Gaussian line profile centered at 
5050 km s~^, the same as observed for shocked Si. No 
absorption with this profile is observed. Figure ^ or ^. 
Since the coUisional ionization rates of Fe II and Si II 
are similar (Lennon et al. 1988), the absence of such 
Fe II absorption implies that there cannot be much 
iron mixed in with the shocked Si. This is consistent 
with the argument in subsection p.5\ that the bulk of 
the shocked Si II must be fairly pure, unmixed with 
other elements. 

On the other hand the observed Fe II profile. Fig- 
ure ^, does suggest the presence of some Fe II mixed 
with unshocked Si, at velocities 7070 km s~^. The 
picture then is that there is Fe mixed with Si at lower 
velocities, but not at higher velocities. This is con- 
sistent with SN la models, such as W7 of Nomoto et 
al. (1984), in which the transition from the inner Fe- 
rich layer to the Si-rich layer is gradual rather than 
abrupt. 

On the blue side of the Fe II profile, if the reverse 
shock is at — 4200kms~^, then shocked Fe II should 
have a Gaussian profile centered at ~ — 2000kms~^, 
with a width comparable to that of the broad red- 
shifted Si features. While some such absorption may 
be present, the similarity between the blueshifted and 
redshifted sides of the Fe II absorption suggests that 
the contribution to the blue side from shocked Fe II 
is not large. This is consistent with expectation from 
the density profile of unshocked Fe II on the red side. 
The mass of Fe II at velocities 4200-7070 kms"^ in- 
ferred from the red side of the profile on the assump- 
tion of spherical symmetry is 0.013 Mq. If this mass 
of Fe II is supposed shocked on the blue side and 
placed at the reverse shock radius of — 4200kms~^, 
the resulting Fe column density is 1.3 x lO^'^cm"^, 
which translates into a peak Fe density of 0.0013 x 
10~^''gmcm~^ at velocity — 2200kms~^, for an as- 
sumed dispersion of 1240 km s~^ like that of the red- 
shifted Si features. This density of shocked Fe is low 
enough that it makes only a minor contribution to the 
Fe II profile in Figure ^. 

In practice, coUisional ionization of shocked Fe II 
reduces its contribution further. For pure iron with 
an initial ionization state of, say, 50% Fe II, 50% 
Fe III (see subsection 3^), we find that the col- 



umn density of shocked Fe II is reduced by a factor 
0.6 to 0.8 X 10"'^^cm~^, which translates into a peak 
Fe II density of 0.0008 x lO^^^gmcm"'^ at velocity 
— 2200kms~^. The shocked column density would 
be even lower if the initial ionization state is higher, 
or if there are other elements mixed in with the Fe, 
since a higher electron to Fe II ratio would make 
coUisional ionization faster. If the initial ionization 
state of the iron is as high as proposed by HF88, 
then the shocked Fe II column density of Fe could 
be as low as 0.1 x lO^^'cm"^, for a peak density of 
0.0001 x 10-24 gmcm-3 a,t velocity -2200 kms'^ in 
Figure |^. 

3.5. Photoionization 

The mass of Fe II inferred here from the Fe II pro- 
file is, according to equation (H), 0.0195±0.0013 Mq 
up to 5000kms-i, and O.O29±O.OO4M0 up to 7070 kms"!. 

Historical and circumstantial evidence suggests that 
SN1006 was a Type la (Minkowski 1966; Schaefer 
1996). Exploded white dwarf models of SNIa predict 
that several tenths of a solar mass of iron ejecta should 
be present, as required to explain SNIa light curves 
(Hoflich & Khokhlov 1996). Thus, as emphasised by 
Fesen et al. (1988) and by HF88, the observed mass 
of Fe II in SN1006 is only a fraction (^ 1/10) of the 
expected total Fe mass. 

Hamilton & Sarazin (1984) pointed out that un- 
shocked SN ejecta will be subject to photoioniza- 
tion by UV starlight and by UV and x-ray emission 
from the reverse shock. Recombination is negligi- 
ble at the low densities here. HF88 presented de- 
tailed calculations of the time-dependent photoioniza- 
tion of unshocked ejecta in SN1006, using deflagrated 
white dwarf models W7 of Nomoto et al. (1984), and 
CDTG7 of Woosley (1987, private communication), 
evolved by hydrodynamic simulation into a uniform 
interstellar medium. HF88 found that in these mod- 
els most of the unshocked Fe was photoionized to 
Fe III, Fe IV, and Fe V. While model W7 produced 
considerably more Fe II than observed in SN1006, 
model CDTG7, which is less centrally concentrated 
than W7, produced an Fe II profile in excellent agree- 
ment with the lUE Fe II 2600 A feature. HF88 con- 
cluded that several tenths of a solar mass of unshocked 
Fe could be present at velocities < 5000 kms^^ in 
SN1006, as predicted by Type la supernova models. 

However, the low ionization state of unshocked Si 
inferred from the present HST observations does not 
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support the high ionization state of Fe advocated by 
HF88. According to the ion fractions given in Table ^, 
unshocked Si is 92 ± 7% Si II. By comparison, HF88 
argued that unshocked Fe is only ~ 10% Fe II. We 
now show that these ionization states of Si II and 
Fe II are not mutually consistent. 

The ionizations of unshocked Si II and Fe II are re- 
lated according to their relative photoionization cross- 
sections, and by the energy distribution of the pho- 
toionizing radiation. Neutral Si and Fe can be ne- 
glected here, since they have ionization potentials be- 
low the Lyman limit, and are quickly ionized by UV 
starlight (Si I in - 20 yr, Fe I in - 100 yr if the 
starlight is comparable to that in the solar neighbor- 
hood), once the ejecta start to become optically thin 
to photoionizing radiation, at ^ 100-200 yr. 

Photoionization cross-sections of Si II and Fe II, 
taken from Reilman & Manson (1979) and adapted to 
include autoionizing photoionization as described by 
HF88, are shown in Figure^. The Figure shows that 
the photoionization cross-section of Fe II is about an 
order of magnitude larger than that of Si II from the 
ionization potential up to the L-shell (autoionizing) 
photoionization threshold of Si II at 100 eV, above 
which the photoionization cross-sections are about 
equal, until the L-shell threshold of Fe II at 700 eV. 
According to HF88 (Table 6, together with Table 2 of 
Hamilton & Sarazin 1984), much of the photoionizing 
emission from the reverse shock is in the UV below 
100 eV. However, there is also some soft x-ray emis- 
sion above 100 eV, which is important for Si II because 
its L-shell photoionization cross-section is larger than 
that of the valence shell. Averaging over the pho- 
toionizing photons tabulated by HF88, we find that 
the effective photoionization cross-section of Fe II is 
about 5 times that of Si II, which is true whether the 
source of the emission in the reverse shock is oxygen, 
silicon, or iron. 

If the effective photoionization cross-section of Fe II 
is 5 times that of Si II, then an unshocked Si II frac- 
tion of 0.92 ± 0.07 would predict an unshocked Fe II 
fraction of (0.92 ± 0.07)^ = OMtoil, considerably 
larger than the desired unshocked Fe II fraction of 
0.1. The 3cr lower limit on the Si II fraction is 0.71, 
which would predict an unshocked Fe II fraction of 
(0.71)^ = 0.18, closer to but still higher than desired. 

A higher ionization state of Fe compared to Si 
might be achieved if the photoionizing emission could 
be concentrated entirely in the UV below 100 eV, 
since then the effective photoionization cross-section 



of Fe II would be about 10 times that of Si II, as il- 
lustrated in Figure ^. This could occur if the pho- 
toionizing emission were mainly from He. In this 
case, the predicted unshocked Fe II fraction would 
be (0.92 ± 0.07)^° = OAStoH, with a 3ct lower limit 
of (0.71)"'^° = 0.03, which is satisfactory. To achieve 
this relatively high level of Fe ionization requires that 
there be little soft x-ray emission from heavy ele- 
ments. It is not clear whether this is possible, given 
the observed x-ray emission from oxygen and silicon 
(Koyama et al. 1995). 

Thus the low ionization state of unshocked Si in- 
ferred in this paper is difficult to reconcile with the 
expected presence of several tenths of a solar mass of 
Fe at velocities <^ 5000 km s^^. Is it possible that the 
unshocked Si is substantially more ionized than we 
have inferred? From Tables and H, the total column 
densities of Si ions, shocked and unshocked, are in the 
ratio 



Nsiu ■■ Nsmi : iVsiiv = 1 : 0.39 : 0.36 



(28) 



If the ionization state of the unshocked Si were this 
high, then a high ionization state of Fe would be al- 
lowed, and the problem would be resolved. In fact, 
in photoionization trials similar to those described by 
HF88, we find that the unshocked Si fractions pre- 
dicted by the CDTG7 model are close to the ratio 



However, in subsection 2.6 we argued both theoret- 
ically and from the observed line profiles that most 
of the observed Si III and Si IV absorption is from 
shocked Si. Might this be wrong, and could in fact 
much or most of the absorption be from unshocked Si? 
And is then our interpretation of the Si II 1260 A pro- 
file as arising mostly from shocked Si II also faulty? If 
so, then much of the tapestry of reasoning in this pa- 
per begins to unravel. For example, we must regard 
as mer ely c oincidental the agreement, found in sub- 
section |2.2| , of the measured parameters of the Si II 
1260 A profile with the energy shock jump condition, 
equations (|1)-(||). We must also conclude that the ob- 
served asymmetry between the red and blucshifted Si 
absorption arises from asymmetry in the initial explo- 
sion, not (or not all) from asymmetry in the ambient 
ISM. For if the blue edge of the redshifted Si II and 
Si IV features (the blue edge of redshifted Si III is 
obscured by Lya) arises from unshocked Si extend- 
ing down to velocities -1-2500 kms~^ (see Figs. |l| and 
^), then there should be, on the assumption of spher- 
ical symmetry, corresponding blueshifted Si absorp- 
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tion outward of —2500 km s ^, which is not seen. 

3.6. Where's the iron? 

Hoflich & Khokhlov's (1996) Table 1 presents a sur- 
vey of 37 SNIa models, encompassing all currently 
discussed explosion scenarios. In their models, the 
ejected mass of ^^Ni, which decays radioactively to 
iron, ranges from 0.10 Mq to 1.07 Mq. Models yield- 
ing 'acceptable' fits to the sample of 26 SNIa consid- 
ered by Hoflich & Khokhlov have ejected ^^Ni masses 
of 0.49-0.83 M0 for normal SNIa, and 0.10-0.18 
for subluminous SNIa. In the subluminous models, 
a comparable amount of Fe is ejected along with the 
5*^Ni (Hoflich, Khokhlov & Wheeler 1995), so the total 
iron ejected in these cases is « 0.2-0.3 Mq. 

In the previous subsection, we argued that the low 
ionization state of unshocked Si inferred from the 
present observations suggests that the ionization state 
of unshocked Fe is also likely to be low. Specifically, if 
the unshocked Si II fraction is Si Il/Si = 0.92 ± 0.07, 
from Table |l|, then the predicted unshocked Fe II frac- 
tion is Fe II/Fe = OMtoH- Correcting the Fe II mass 
of Mpcii = 0.029±0.004Mq up to 7070kms-\ equa- 
tion (p6|), for the ionization state of the Fe yields a 
total inferred Fe mass of A/pe = 0.044lo;o^3 Mq up to 
7070 km s-\ with a 3a upper limit ofMpo < 0.16Mq. 
These Fe masses are lower than predicted by either 
normal or subluminous models of SN la. 

A low ionization state of Fe is supported by the 
HUT observations of Blair et al. (1996), who looked 
for Fe III 1123 A absorption in the background SM 
star. If there is a large mass of Fe in SN1006, signifi- 
cantly larger than the observed Fe II mass, then cer- 
tainly Fe III should be more abundant than Fe II, and 
from detailed models HF88 predicted Fe III/Fe II = 
2.6. Blair et al.'s best fit is Fe III/Fe II = l.I ± 0.9, 
and their 3f7 upper limit is Fe III/Fe II < 3.8. This 
result does not support, though it does not yet defi- 
nitely exclude, HF88's prediction. 

Neither of the above two observational evidences 
favoring a low ionization state of Fe, hence a low mass 
of Fe in SN1006, is yet definitive. To settle the issue 
will require re-observation of the Fe III 1123 A line 
at a higher signal to noise ratio. The Far Ultraviolet 
Space Explorer (FUSE) should accomplish this. 

4. Worries 

In this paper we have attempted to present a con- 
sistent theoretical interpretation of the broad Si and 



Fe absorption features in SN1006. While the overall 
picture appears to fit together nicely, the pieces of the 
jigsaw do not fit perfectly everywhere. In this Section 
we highlight the ill fits. 

What causes the discrepancy between the profiles 
of the redshifted Si II 1260 A and Si II 1527 A fea- 
tures? This discrepancy was originally pointed out 
and discussed for the lUE data by Fesen et al. (1988), 
and the discrepancy remains in the HST data (WCH- 
FLS96). The discrepancy is especially worrying for 
the present paper because the excess in the Si II 
1260 A profile compared to Si II 1527 A (see WCH- 
FLS96, Figure 2) looks a lot like what we have inter- 
preted here as the unshocked component of the Si II 
1260 A absorption (Figure g). 

We have argued that the redshifted Si II, Si III, and 
Si IV absorption is caused mostly by shocked Si, yet 
it is not clear that the observed relative column den- 
sities are naturally reproduced in coUisionally ionized 
gas. Specifically, one might expect relatively more 
Si III, or relatively less Si II or Si IV in the shocked 
Si. On the other hand the observed relative column 
densities of Si are naturally reproduced in unshocked, 
photoionized gas. Is our interpretation at fault? 

The best fit dispersion of the redshifted Si IV 
feature. Figure ^, is 1700 ± lOOkms"^, which is 
4.5(7 larger than that of the Si II 1260 A feature, 
1240 ± 40kms^^. What causes this discrepancy? 

Docs the d ensi ty of shocked Si vary a little or a lot? 
In subsection 2.1 we argued that the Gaussian profile 
of the shocked Si suggests little temperature varia- 
tion, hence little density variation. On the other hand 
the unshocked Si density profile below 7070 km s~^ 
(Fig. |l|) suggests a density profile increasing steeply 
outwards, and there are other clues hinting at the 
same thing: the large column density of shocked Si II, 
subsection 2.5, and the need for a high density to ion- 
ize Si on the near side above Si IV, subsection 2.9. Is 
there an inconsistency here? 



In subsection 2.5 we argued that the high observed 
column density of shocked Si II indicates a low mean 
electron to Si II ratio, ng/'^-siii ^ 1-3, equation (p^. 
Higher ratios would cause more rapid coUisional ion- 
ization of Si II, reducing the column density below 
what is observed. This limit on the electron to Si II 
ratio is satisfactory as it stands (the limit is some- 
what soft, and at least it is not less than 1), but it is 
uncomfortably low, making it difficult to admit even 
modest quantities of other elements, such as sulfur. 
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which might be expected to be mixed with the siHcon. 
Is there a problem here, and if so have we perhaps 
overestimated the contribution of shocked compared 
to unshocked Si in the Si II absorption? 



In subsection 2.9 we estimated that if the reverse 
shock on the near side of SN1006 is at — 4200kms~^, 
then under a 'simplest' set of assumptions there 
should be an observable column density of blueshifted 
Si IV, contrary to observation. However, we also 
showed that the predicted column density is sensi- 
tive to the assumptions, and that it is not difficult to 
bring the column density of Si IV below observable 
levels. Is this explanation adequate, or does the ab- 
sence of blueshifted Si absorption hint at asymmetry 
in the supernova explosion? 

Is the sharp blue edge on the Fe II features at 
— 4200kms~^ real? Further observations of the Fe II 
features at higher resolution would be helpful in de- 
ciding this issue. 

Notwithstanding our reanalysis of the Fe II fea- 
tures, there remains some suggestion of high velocity 
blueshifted absorption outside — 4200kms~^, perhaps 
to — 7000kms~^ or even farther. Figures || and |[ Is 
this absorption real? If so, then the arguments of 
subsection 2.8 fail, and the absence of blueshifted Si 



absorption must be attributed to intrinsic asymmetry 
in the initial supernova explosion. 

F inal ly, there is the problem discussed in subsec- 
tion 3.6: where is the iron? 



5. Summary 

We have presented a consistent interpretation of 
the broad Si and Fe absorption features observed in 
SN1006 against the background SM star (Schweizer 
& Middleditch 1980). 

We have argued that the strong redshifted Si II 
1260 A absorption feature arises from both unshocked 
and shocked Si, with the sharp red edge of the feature 
at 7070 km s""'^ representing the free expansion radius 
of the reverse shock on the far side of SN1006, and 
the Gaussian blue edge signifying shocked Si (Fig. |^). 
Fitting to the Si II 1260 A line profile yields three 
velocities, the position of the reverse shock, and the 
velocity and dispersion of the shocked gas, permit- 
ting a test of the energy jump condition for a strong 
shock. The measured velocities satisfy the condition 
remarkably well, equations (I)-®. The Si II 1260 A 
line thus provides direct evidence for the existence of 
a strong shock under highly coUisionless conditions. 



The energy jump condition is satisfied provided 
that virtually all the shock energy goes into ions. This 
evidence suggests little or no coUisionless heating of 
electrons in the shock, in agreement with recent evi- 
dence from UV line widths and strengths (Raymond 
et al. 1995; Laming et al. 1996). 

The observed column density of shocked Si II is 
close to the column density expected for steady state 
collisional ionization behind a shock, provided that 
the electron to Si II ion ratio is low. From the low elec- 
tron to Si II ratio, we have argued that the shocked 
Si is probably of a fairly high degree of purity, with 
little admixture of other elements. More directly, the 
absence of Fe II absorption with the same line pro- 
file as the shocked Si indicates that there is little Fe 
mixed with the shocked Si. On the other hand, there 
is some indication of absorption by Fe II at the ve- 
locity 5600-7070 kms~^ of the unshocked Si II, which 
suggests that some Fe is mixed with Si in the lower 
velocity region of the Si layer. 

We have proposed that the ambient interstellar 
density on the far side of SN1006 is anomalously 
low compared to the density around the rest of the 
remnant, so that the remnant bulges out on the far 
side (Fig. 0). This would explain several observa- 
tional puzzles. Firstly, it would explain the absence of 
blueshifted Si absorption matching the observed red- 
shifted Si absorption. If the interstellar density on the 
near side is substantially larger than on the far side, 
then the reverse shock on the near side would be fur- 
ther in, so that all the Si on the near side could have 
been shocked and coUisionally ionized above Si IV, 
making it unobservablc in absorption. Secondly, if the 
velocity on the far side is anomalously high because of 
the low interstellar density there, it would resolve the 
problem noted by Wu et al. (1983) and subsequent 
authors of how to fit gas expanding at ~ 7000 km s~^ 
within the confines of the interstellar shock. Thirdly, 
a low density on the far side would explain how the 
5050 kms~^ velocity of shocked Si there could be so 
much higher than the 1800 kms^^ velocity (assuming 
no coUisionless electron heating) of gas behind the in- 
terstellar shock along the NW filament (Smith et al. 
1991; Raymond et al. 1995). Finally, the density of Si 
on the far side inferred from the Si absorption profiles 
is one or two orders of magnitude too low to yield Si 
x-ray emission at the observed level (Koyama et al. 
1995). Again, an anomalously low density on the far 
side is indicated. 

The notion that the reverse shock on the near side 
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has moved inward much farther, to lower velocities, 
than on the far side conflicts with our earlier con- 
clusion (WCFHS93) that there is blueshifted Fe II 
absorption to velocities ~ — 8000 km s^^. Reanalyz- 
ing the Fe II data, we find that the evidence for such 
high velocity blueshifted Fe II absorption is not com- 
pelling. In the WCFHS93 analysis, the main evidence 
for high velocity blueshifted Fe II comes from the Fe II 
2383, 2344, 2374 A feature. However, the 2344 A com- 
ponent, which lies at — 4900kms~^ relative to the 
principal 2383 A component, confuses interpretation 
of the blue wing of the feature. The Fe II 2600, 2587 A 
feature is cleaner, and it shows a sharp blue edge at 
—4200 kms""'^, which we interpret as representing the 
free expansion radius of the reverse shock on the near 
side of SN1006. In our reanalysis of the Fe II fea- 
tures, we adopt a rigorous approach to the subtrac- 
tion of narrow interstellar and stellar lines, requiring 
that lines subtracted have the correct positions and 
dispersions, and have mutually consistent strengths. 
In particular, we subtract the narrow Fe II 2344 A line 
with a strength consistent with the other narrow Fe II 
lines, which strength is substantially greater than the 
apparent strength. Subtraction of this line introduces 
a sharp blue edge on the deconvolved Fe II 2383, 2344, 
2374A feature at the same place, « — 4200kms~^, as 
the Fe II 2600, 2587 A feature. The resulting decon- 
volved Fe II profiles (Fig. |^) are in good agreement 
with each other. 

The mass and velocity distribution of Si and Fe 
inferred in this paper provides useful information for 
modeling the remnant of SN1006 (see Fig. |^). Freely 
expanding unshocked Si on the far side extends from a 
low velocity of 5600 ± 100 kms^^ up to the position of 
the reverse shock at 7070 kms^^. Above this velocity 
the Si is shocked, and information about its detailed 
velocity distribution before being shocked is lost. The 
total mass of Si, both unshocked and shocked, inferred 
from the Si II, Si III, and Si IV fines is Msi = 0.25 ± 
0.01 Mq, on the assumption of spherical symmetry. 

We have argued that the observed broad Fe II ab- 
sorption arises almost entirely from unshocked freely 
expanding Fe. The mass of Fe II inferred from 
the cleaner, red side of the mean Fe II profile is 
Mpeii = 0.0195 ± 0.0013 M0 up to 5000 kms^^, and 
Mpeii = 0.029±0.004Mq up to 7070 kms^^, again on 
the assumption of spherical symmetry. These masses 
include a small positive adjustment (0.003 Mq) re- 
sulting from the offset of the SM star from the pro- 
jected center of the remnant. 



Our analysis of the Si lines indicates a low ion- 
ization state for the unshocked silicon, with Si Il/Si 
= 0.92 ± 0.07. Such a low state would imply a 
correspondingly low ionization state of unshocked 
iron, with Fe II/Fe = 0.66;o 22- If this is cor- 
rect, then the total mass of Fe up to 7070 km s~^ 
is Afpc = 0.044tQ;oi3 Mq with a 3cr upper limit of 
A/pc < 0.16 Mq. The absence of Fe II absorption 
with a profile like that of the shocked Si II suggests 
that there is not much more Fe at higher velocities. 
Such a low mass of Fe conflicts with the expectation 
that there should be several tenths of a solar mass of 
Fe in this suspected Type la remnant. A low ioniza- 
tion state of Fe and a correspondingly low Fe mass is 
consistent with the low Fe III/Fe II — 1.1 ± 0.9 ratio 
measured by Blair et al. (1996) from HUT observa- 
tions of the Fe III 1123 A line in the spectrum of the 
SM star. However, neither the present observations 
nor the HUT data are yet conclusive. 

Re-observation of the Fe III 1123 A line at higher 
signal to noise ratio with FUSE will be important in 
determining the ionization state of unshocked Fe in 
SN1006, and in resolving the question, Where's the 
iron? 
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respondence on the HUT data, and Graham Parker 
and Mike ShuU for advice on respectively stellar and 
interstellar lines. Support for this work was provided 
by NASA through grant number GO-3621 from the 
Space Telescope Science Institute, which is operated 
by AURA, Inc., under NASA contract NAS 5-26555. 



A. Ionization times in non-steady state 

This Appendix demonstrates the assertion made 



in subsection 2.5, that the column density per in- 



terval of ionization time is reduced compared to the 
steady state value if the density of shocked gas in- 
creases dow nstr eam, and vice versa. The main result 
is equation (A7). 

Define a collisional ionization time t (not to be con- 
fused with optical depth) by the integral over time 
of the number density n in a Lagrangian gas ele- 
ment, from the time when the gas element was 
first shocked, to the present time t: 



idt 



(Al) 
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The density n in the definition (Al) of the ionization 
time is the ion density. One might think to use the 
electron density Ue rather than the ion density, since 
the coUisional ionization is by electron impact, but 
in fact the electron to ion ratio is itself determined 
by the ionization time, for any given initial ionization 
state, so there is no loss of generality in the definition 
®. 

It is useful to regard the shock age ts as a La- 
grangian coordinate labeling different gas elements. 
Then the ionization time r of different gas elements 
at fixed time t varies according to 



dr 




/■* dn 




t 


Ju dtl 



dt 



(A2) 



where Ug is the postshock density at the time the gas 
was shocked. The minus sign in front of n^. reflects 
the fact that gas elements farther downstream were 
shocked at earlier times. If the pressure is spatially 
constant, which should be a good approximation, and 
if the gas density varies adiabatically as the pressure 
changes in time, then the densities n in different el- 
ements vary with time in proportion to each other. 
It follows that the logarithmic derivative of density 
n with respect to the Lagrangian coordinate ts is in- 
dependent of time. It is convenient to denote this 
dimensionless logarithmic derivative by 7: 



l{ts) 



d\nn 



(?lnt» 



(A3) 



which is positive if the density at fixed time t increases 
downstream. Then equation (A2) becomes 



dr 



9 In ts 



7T 



(A4) 



Now consider the column density N contributed by 
gas elements shocked at different times ts- The rate 
at which mass M of ions mass m enters the shock at 
time i,, when the shock radius is r„ 



M = 



dM 
dts 



IS 



= 47rmr^nP™''''ws = ■nmr'lnsVs (A5) 



the postshock density being 4 times the preshock den- 
sity for a strong shock, Ug = 4nP''°^'^. The minus sign 
—dM/dtg here arises because we choose to label the 
interior mass M(ts) so that it increases downstream 
of the shock, in the direction of decreasing ts. Mass 
is conserved in Lagrangian elements, and the column 



density N at time t when the shock radius is r varies 
with t^ as 



dN 



d\nt^ 



1 



dM 



ri UsVsts 



Aumr'^ dint. 



(A6) 



which is diluted by the geometric factor r^/r^ com- 
pared to the column density at the time ts when the 
gas was shocked. 



Comparing the column density (A6) to the ion- 
ization time (A4) shows that the column density per 
interval of ionization time observed at a time t is 



dN 
'd7 



1 



4r2 [l+7r/(nA)] 



(A7) 



which is the principal result of this Appendix. In 



steady state, equation (A7) reduces to dN/dT\t = 
Vs/i. Compared to steady state, the non-steady state 
column density per unit ionization time (A7) thus 
contains two factors discussed in subsection 2.5, the 



geometric dilution factor rj jr , and the density pro- 
file factor [1 -|- 7r/(nsts)]~^. Since r and Usts are al- 
ways positive, this latter factor is less than or greater 
than one depending on whether 7, equation ( |A3| ), is 
greater than or less than zero. In other words, if the 
density n increases downstream, 7 > 0, then the col- 
umn density is reduced compared to steady state, and 
vice versa, as was to be proved. 
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Fig. 1.— HST G130H spectrum relative to the stel- 
lar continuum around the redshifted Si II 1260.4221 A 
feature, showing the best fit Gaussian profile of 
shocked Si II, and the residual unshocked Si II. Up- 
per axis shows velocity in the rest frame of the Si II 
line. Measured parameters of the feature are given 
in Table |^. As elsewhere in this paper, we assume a 
stellar continuum which is linear in log F-\og A, and 
fit the continuum to ostensibly uncontaminated re- 
gions around the line. The uncertainties in the pa- 
rameters given in Table ^ include uncertainty from 
placement of the continuum. The adopted stellar con- 
tinuum is logF ~ log(4.1 X 10^^'* ergs~^ cm~^ A ) — 
2.31og{A/1260A). 
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Fig. 2.— HST G130H spectrum (sohd hue) rela- 
tive to the adopted stellar continuum around the red- 
shifted Si III 1206. 500A (/ = 1.669) feature. Upper 
axis shows velocity in the rest frame of the Si III line. 
The position and width of the fitted Gaussian profile 
of Si III (dotted line) have been constrained to be the 
same as that of the Si II 1260 A feature. Also shown 
is the residual (dashed line) after subtraction both of 
the Si III fitted Gaussian, and of the contribution of 
redshifted Si II 1193, 1190 A (dotted line) assumed to 
have the same profile as the Si II 1260 A feature. The 
residual shows, besides Ly a emission and absorption, 
a weak indication of absorption by unshocked Si III 
at velocities 5500-7000 kms~^. The adopted stellar 
continuum is the same as that for the Si II 1260 A 
feature in Figure El 
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Fig. 3.— HST G130H spectrum relative to the 
adopted stellar continuum around the redshifted Si IV 
1393.755, 1402. 770A (/ = 0.5140, 0.2553) feature. 
Upper axis shows velocity in the rest frame of the 
weighted mean wavelength of the feature. Dotted 
lines show the best fit to a Gaussian pair with line 
center and dispersion constrained to be that of Si II 
1260 A for each component of the doublet, and with 
optical depths fixed equal to the ratio of oscillator 
strengths times wavelengths of the doublet. The best 
fit dispersion of the Si IV is 1700 ± 100 kms^^, which 
is 4.5(7 larger than the 1240 kms^^ dispersion of the 
fit shown here. There is no evidence of unshocked 
Si IV. The adopted stellar continuum is logi^ = 
log(2.9xl0-i''ergs-icm-2 A"^)-1.41og(A/1397A). 



Fig. 6. — Inferred density profile of ejecta in SN1006. 
The Fe II profile is from the mean of the two decon- 
volved broad Fe II absorption features shown in Fig- 
ure |. The dotted line at +5600-7070 kms^^ is the 
profile of unshocked Si II from the redshifted Si II 
1260 A absorption feature in Figure ^ The dashed 
line above 7070 km s^^ is a plausible extrapolation 
of the total Si density before it was shocked: it is 
a quadratic function of velocity, equation (|2^), which 
approximately reproduces the observed profile of un- 
shocked Si II, and which contains the observed total 
Si mass of 0.25 Mq (assuming spherical symmetry), 
equation (p7|). 
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Fig. 7. — Schematic diagram, approximately to scale, 
of the structure of the remnant of SN1006 inferred in 
this paper. The remnant on the far side bulges out 
because of the low interstellar density there. Shaded 
regions represent silicon ejecta, both shocked and un- 
shocked. Iron, both shocked and unshocked, lies in- 
side the silicon. The background SM star is offset 
slightly from the projected center of the remnant. 



00 



8 



o 



a 

N 

'S 
o 



o 



.01 



I I 1 1 m 




10 100 1000 

Photon energy (eV) 



10000 



Fig. 8. — Photoionization cross-sections of Si II 
and Fe II (Reilman & Manson 1979, adapted to in- 
clude autoionizing photoionization as described by 
HF88). The ionization state of unshocked ejecta de- 
pends on photoionization, and the plotted photoion- 
ization cross-sections are important in relating the 
ionization state of unshocked iron, characterized by 
the ratio Fe II/Fe, to the ionization state of unshocked 
silicon, characterized by the ratio Si 11/ Si. 
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